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ABSTRACT 

The photospheres of about 10-20% of main sequence A- and B-type stars exhibit a 
O \ wide range of chemieal pecuharities, often associated with the presence of a magnetic 

field. It is not exactly known at which stage of stellar evolution these chemical pe- 
culiarities develop. To investigate this issue, in this paper we study the photospheric 
compositions of a sample of Herbig Ae and Be stars, which are considered to be the 
pre-main sequence progenitors of A and B stars. We have performed a detailed abun- 
dance analysis of 20 Herbig stars (three of which have confirmed magnetic fields), 
^ . and one dusty young star. We have found that half the stars in our sample show A 

Boo chemical peculiarities to varying degrees, only one star shows weak Ap/Bp pe- 
culiarities, and all the remaining stars are chemically normal. The incidence of A Boo 
chemical peculiarities we find in Herbig stars is much higher than what is seen on the 
main sequence. We argue that a selective accretion model for A Boo star formation is 
a natural explanation for the remarkably large number of A Boo stars in our sample. 
We also find that the magnetic Herbig stars do not exhibit a range of chemical compo- 
sitions remarkably different from the normal stars: one magnetic star displays A Boo 
chemical peculiarities (HD 101412), one displays weak Ap/Bp peculiarities (V380 Ori 
A), and one (HD 190073) is chemically normal. This is completely different from what 
, is seen on the main sequence, where all magnetic A and cool B stars show Ap/Bp 

. chemical peculiarities, and this is consistent with the idea that the magnetic field 

■ precedes the formation of the chemical peculiarities typical of Ap and Bp stars. 

Key words: stars: magnetic fields, stars: abundances, stars: chemically peculiar, 
stars: pre-main-sequence 
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1 INTRODUCTION 

In main sequence A and B stars a wide range of peculiar 
chemical abundances are observed. While much study has 
been devoted to these peculiar objects, a lot remains to be 
learnt about the physical processes behind these chemical 
peculiarities. An avenue of study that has not yet been prop- 
erly investigated is to examine the chemistry in the pre-main 
sequence progenitors of main sequence peculiar stars. 

Herbig Ae and Be (HAeBe) stars are pre-main se- 
quence A and B stars, and thus evolve into the wide 
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range of different main sequence A and B stars. Obser- 
vationally, Herbig stars are generally identified by optical 
emission lines, infrared excess, and are usually associated 
with nebulous regions (.Vieira et al.. 2003i ). Chemically, these 
stars are thought to posse ss approximately solar abundances 
(lAcke fc Waelkeng |2004| ). similar to most young, nearby 
main sequence stars. 

A small number of HAeBe stars have recently been 
fo und to have str o ng, globally orde r ed magnetic fields 



i Donati et all 1199 



._ , -..^^l f I^OO'^ ' ICatala et all 

l2007l : lAlecian et all iToOSal M . The strength and morphoh 
ogy of these magnetic fields are very similar to those 
seen in the magnetic chemically peculiar Ap and Bp 
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stars (|Alecian et al.ll200 8a': 'Folso m et al]|2008l : lAlecian et alj 
|2009| ). making these magnetic HAeBe stars st rong candi- 
dates for the progenitors of Ap and Bp stars (|Wade et alJ 
This discovery in particular raises the question of 
whether any sign of chemical peculiarity can be found on 
the pre-main sequence. 

In main sequence A and cooler B stars, magnetic fields 
are always seen together with t he characteristic chem ical pe- 
culiarities of Ap and Bp stars (jAuriere et al.|[2007^ . Hotter 
magnetic B stars also usually show chemical peculiarities, 
such as He-weak or He-strong stars. Consequently, if mag- 
netic HAeBe stars evolve into Ap and Bp stars, then at 
some point they must develop chemical peculiarities. De- 
tecting such peculiarities on the pre-main sequence would 
provide new constraints on the timescales on which, and the 
conditions under which, these peculiarities evolve. 

A few detections of chemical peculiarities in very young 
stars near or on the zero-age main sequence provide tanta- 
lising hints that chemical peculiarities may be common in 
magnetic Herbig Ae/Be stars. Some notable ca ses of this are 



V380 Ori (|Alecian et al.1 



20091), HP 72106 A jFolsom et all 



I2008D . NGC 661 1 W601 (lAlecian et all2008bD . and perhaps 
NGC 2244-334 (|Bagnulo et al.ll2004f ) whi ch is a very young 



main sequence star. Recent modelling bv lVick et all |201li 



of chemical diffusion in the presence of modest mass loss 
suggest it is possible for chemical peculiaritie s to form dur- 
ing th e pre-main sequence phase. Interestinglv. lCowlev et al] 
l|2010t ) found a different form of chemical peculiarities, A 
Bootis peculiarities, in the Herbig Ae star HD 101412. 

A Bootis stars are mostly main sequence A stars, with 
strong underabundances of many metals, particularly iron 
peak elements. Lighter elements, specifically C, N, O, S, 
and in some cases Na, have normal abundances, while 
iron peak elements are usually depleted by ^1 dex (e.j 



IVenn fc Lambertll990l : lAndrievskv et al.ll2002l : lHeiteil200: 

This is in contrast to almost all other chemically peculiar 
A and B stars, which are characterised by strong over- 
abundances of iron peak elements. A Boo stars have not 
been found to have magnetic fi elds, unlike Ap and Bp stars 
l|Bohlender fc Landstreetlll990l ). A Boo stars have a distri- 
bution of V s in i values that is the same as for normal A and 
B stars (e.g. lAbt fc Morrelllll995! ). unlike Ap or Am stars. 

The cause of the peculiar abundances seen in A Boo 
stars remains unknown, though a number of hypotheses 
have been suggested. The most popular hypothesis is that 
the peculiarities are a result of a selective ac cretion pro- 
cess llVenn fc Lambertll 19901 : IWaters et af]|l992l ') . In this sce- 
nario, gas depleted in heavier elements is accreted prefer- 
entially, building up a layer of relative underabundance in 
heavier elements at the surface of the star. A proposed mech- 
anism for this suggests that heavier elements are bound into 
dust grains, which are then preferentially driven away from 
the star by its radiation, while gas, which is depleted in heav- 
ier elements as a con sequence of dust formation, is more 
readily accreted (e.g. lAndrievs kv fc PaunzenI bOOOj ). Thus, 
because of the lack of significant convective mixing in the 
atmospheres of A-type stars, a layer that is underabundant 
in heavier elements is quickly built up at the surface of the 
star. 

A variation on this hypothesis is that the A Boo pecu- 
liarities result from selective accretion of gas in the interstel- 
lar medium, rather than pre-existing circumstellar material 



iKamp fc Paunzenl l2002''). In this scenario the star passes 
through a diffuse interstellar cloud, in which heavier ele- 
ments are already preferentially bound into dust grains. The 
star then accretes the metal poor gas, while driving away the 
dust by radiation pressure. 

Other hypotheses have been put forward to explore 
the origins of A Boo stars. For example it has been sug- 
gested that A Boo stars are actua lly spectroscopic bi- 
naries (|Faraggiana fc Bonifacio! |2005| ). thus the apparent 
under abundances would be due to continuum emission 
from a secondary, but t his hypothesis is largely discounted 
(|Stiitz fc PaunzenI I2OO6I ). Radiatively driven atomic diffu- 
sion, which is the most widely accepted mechanism for form- 
ing peculiarities in most chemically peculiar A and B stars, 
has been c onsidered for A Boo stars with the addition of 
mass loss (|Michaud fc Charlandl [T986I ) . However, any tur- 
bulent or rotational mixing impedes the proc ess sufficiently 
to pr event it from forming A Boo peculiarities (|Charbonneaul 
119931 ). Since A Boo stars rotate with similar v sin i to nor- 
mal A stars, they potentially can have significant rotational 
mixing. Thus atomic diffusion is generally considered to be 
insufficient to cause A Boo peculiarities. 

Numerical models of A Boo star atmospheres includ- 
ing the accretion of material depleted in heavier ele- 
me nts, together with atomi c diffu sion, have bee n mad e 
by iTurcotte fc Charbonneeail (|l993l ) and iTurcottd l|2002h . 
They find that those models can produce A Boo peculiar- 
ities quickly (~0.1 Myr) for gas accretion rates of ^lO"'^^ 
M0yr~^, but the peculiarities in these models do not last 
long after accretion has stopped (~1 Myr). The peculiarities 
in these models appear to persi st despite rotational mixing 
(|Turcotte fc Charbonneaulll993l ). but larger surface convec- 
tion zones re quire larger accr etion rates to produce the same 
pecuharities (|Turcottell200^ . 

In order to help address the question of A Boo star for- 
mation, and to investigate when Ap/Bp stars develop chem- 
ical peculiarities, we have analysed high resolution spectra 
of 20 Herbig Ae/Be stars, and one dusty young star. We 
have determined atmospheric parameters and performed a 
detailed abundance analysis for these stars, with the goal of 
detecting chemical peculiarities. Three of these stars have 
confirmed magnetic field detections, and thus may be the 
progenitors of Ap/Bp stars. The other 18 stars have been 
carefully searched for magnetic fields, but no confirmed mag- 
netic fields have been found. 



2 OBSERVATIONS 

Observations for this study were obtained with the Echelle 
Spectropolarimetric Device for the Observation of Stars (ES- 
PaDOnS) instrument at the Canada France Hawaii Tele- 
scope (CFHT). ESPaDOnS is a high-resolution echelle spec- 
tropolarimeter, providing nearly continuous wavelength cov- 
erage from 3700 to 10500 A at a resolution of 7? = 65000. 
All observations were obtained in spectropolarimetric mode, 
which provides circularly polarised Stokes V spectra, as well 
as total intensity Stokes I spectr a. Data reduction w as per- 
formed with the Libre-ESpRIT (|Donati et al.lll997l ) pack- 
age, which is optimised for the ESPaDOnS instrument, and 
performs calibrations and optimal spectrum extraction. 

The observations were obtained over several years, be- 
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tween 2004 and 2006 as part of an extended campaign 
investigating the presence of magnetic fields in Herbig 
Ae and Be s t ars. Tl iis study is discussed as a whole by 
lAlecian et al. ( 2012a ). and individua l resul ts are reported 



by 



Wade et al 



(l2005h. ICatala et all (120071'). lAlecian et all 



ll2008al). lAlecian et all (|2008bl ). iFolsom et al.l (|2008l ). and 



lAlecian et al.l ( 20091 ). The observations that were analysed 
for this pap er are reported in Ta ble [T] Additional obser- 
vations from lAlecian et all (|2012al ) were used to check for 
variability in spectral lines. The high resolution and high 
S/N necessary for the detection of magnetic fields in metal- 
lic lines also makes these observations very well suited to 
a detailed abundance analysis, thus we have an excellent 
dataset for our study. 

The stars analysed in this s tudy are only a subset of the 
full set of targets observed by lAlecian et al.l (|2012al ). The 
stars selected for our study were chosen to cover a range 
of Tefi and v sin i values. The stars were also selected to 
generally have only modest amounts of emission in their op- 
tical spectra, which allows for a more complete and accurate 
abundance analysis. Preference was given to stars that were 
not clearly spectroscopic binaries, and to observations with 
a peak S/N over 200. Stars with a confirmed magnetic field 
and a Tcff below 15000 K were all incl uded (except for HP 
72106 A, which was studied in detail bv lFolsom et al.ll2008l ). 

Observations of one star (HD 101412) were obtained 
from the Anglo-Australian Telescope (AAT) with the Uni- 
versity College London Echelle Spectrograph (UCLES) to- 
gether with SEMPOL polarimeter. This instrument consists 
of a bench mounted cross-dispersed echelle spectrograph 
(UCLES) fibre fed from a Cassegrain mounted polarimeter 
unit (SEMPOL), and is fundamentally a similar instrument 
to ESPaDOnS. De tails of the instrument can be found in 
ISemel et all (|l993l ) and iDonati et all (|l99ll2003D . and the 
observing run that the o bservations of HD 10 1412 were ob- 
tained in is described bv lMarsden et aP (j201ll ). Data reduc- 
tion and optimal spectral extr action were performed with a 
generic version of ESpRIT (i Donati et al.| [l997). 

Supplemental archival observations from the FORSl in- 
strument at the Very Large Telescope (VLT) were used for 
examining Balmer lines in many stars, and are listed in Ta- 
ble (2] FORSl is a low resolution spectropolarimeter, which 
produces an observed spectrum entirely in one order. This 
is useful as it reduces the possibility of normalisation errors 
across Balm er lines. The obser vations were obtained from 
the study of I Wade et all (|2007l ). and a detailed description 
of the observations and data reduction techniques can be 
found therein. While IWade et all (|2007l ) focused on detect- 
ing magnetic fields, we only concern ourselves with Balmer 
line profile shapes. 



Table 1. Table of observations for which an abundance analy- 
sis was performed. The Heliocentric Julian date, total integra- 
tion time, and peak signal-to-noise ratio (per 1.8 kms~^ spectral 
pixel) in Stokes I are presented for each observation. 



Object 



HJD 



Integration Peak S/N 









Time (s) 


I 


HD 17081 


2453422 


7286 


480 


878 


HD 31293 


2453423 


8512 


2400 


471 


HD 31648 


2453423 


8840 


2400 


334 


HD 36112 


2453421 


7718 


2400 


239 


HD 68695 


2453423 


9585 


2400 


112 


HD 139614 


2453422 


0790 


3600 


222 


HD 141569 


2454167 


0581 


5400 


877 


HD 142666 


2453424 


0544 


3600 


288 


HD 144432 


2453423 


1169 


3200 


312 


HD 163296 


2453607 


7370 


1200 


538 


HD 169142 


2453606 


8006 


2000 


394 


HD 176386 


2453607 


7587 


1600 


479 


HD 179218 


2453608 


8447 


1600 


532 


HD 244604 


2453607 


0733 


3600 


276 


HD 245185 


2453421 


8207 


4800 


146 


HD 278937 


2453422 


7664 


4800 


167 


T Ori 


2453607 


1183 


3600 


204 


HD 101412 


2454194 


9520 


7200 


150 


HD 190073 


2454167 


1543 


2700 


496 


V380 Ori 


2453609 


0940 


4800 


256 



Table 2. FORSl observations used for Balmer line fitting, ob- 
tained by Wade et al. (2007). The Heliocentric Julian date, total 
integration time, and peak signal-to- noise ratio (per 64 kms~^ 
spectral pixel) in Stokes / are presented for each observation. 



Object 



HJD 



Integration Peak S/N 









Time (s) 


/ 


HD 17081 


2453331 


5530 


12 


3800 


HD 31293 


2453331 


6774 


140 


4050 


HD 31648 


2453331 


6929 


319 


4225 


HD 36112 


2453331 


7097 


320 


3725 


HD 68695 


2453332 


7761 


640 


1950 


HD 141569 


2453062 


8424 


1295 


5650 


HD 142666 


2453063 


8535 


2800 


2775 


HD 144432 


2453062 


8983 


2100 


4575 


HD 244604 


2453331 


6217 


1600 


4150 


HD 245185 


2453331 


6527 


1600 


3925 


HD 278937 


2453330 


6795 


1860 


2550 


T Ori 


2453332 


6371 


1440 


1900 


HD 101412 


2453062 


7977 


1350 


1575 


HD 190073 


2453330 


5153 


215 


3700 


V380 Ori 


2453330 


7576 


480 


1450 



3 FUNDAMENTAL PARAMETERS 

Temperature and surface gravity were first measured by fit- 
ting Balmer lines, and then confirmed and, if possible, re- 
fined by enforcing local thermodynamic equilibrium (LTE) 
ionisation and excitation balances in the metallic line spec- 
tra of the stars during the modelling procedure. 

Balmer line fitting was performed on FORSl spectra 
when available, since it is easier to normalise broad spec- 
tral features when they are recorded entirely in one spec- 



tral order. When FORSl spectra were not available, ES- 
PaDOnS spectra were used and careful attention was paid 
to the normalisation of Balmer lines. Continuum points well 
outside the wings on both sides of a Balmer line were care- 
fully chosen, a low order polynomial was then fit through 
these points, and the polynomial was checked to ensure it 
varied only a small amount over the Balmer line. In the 
ESPaDOnS spectra, careful data reduction and the absence 
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of instrumental artifacts produced an excellent correspon- 
dence between spectral orders in regions where the orders 
overlap. These overlap regions were used to evaluate and 
minimise errors in normalisation across spectral orders. For 
stars with both FORSl and ESPaDOnS spectra available 
the normalised Balmer lines were compared to ensure the ac- 
curacy of the normalisation procedure for ESPaDOnS spec- 
tra. A good agreement between the normalised ESPaDOnS 
and FORSl Balmer lines was consistently found. In all cases 
potential normalisation errors are included in the quoted un- 
certainties of physical parameters. 

Careful attention was paid to avoid contamination of 
the Balmer lines by emission. In our observations the cores of 
most Balmer lines were partially infilled by emission. Ha was 
often heavily contaminated by emission, and in many stars 
was avoided entirely. Despite this, the wings of H7 and 
generally remain free from emission. The wavelength range 
of possible contamination in the higher Balmer lines can usu- 
ally be assessed by the more prominent emission in Ho. As 
a consequence of these complications, the fitting procedure 
was performed by hand and rather conservative uncertain- 
ties were adop ted. Synthetic Balmer lines were computed 
from ATLAS9 l|Kurucj|l993t ) model atmospheres using so- 
lar chemical abundances. 

Problematically, for many of our stars there is a sub- 
stantial degeneracy in effective temperature and surface 
gravity when fitting only the wings of Balmer lines. As a 
result there is a large covariance between T^s and log g, 
and our uncertainties often represent long ellipses in the 
two dimensional parameter space. Typically these degener- 
ate regions have a slope of 4-250 K in Tcff for -f 0.1 dex in 
log g, though the slope generally becomes steeper towards 
hotter temperatures. 

When good quality fits to the observations were obtain- 
able, fitting of Teff and log g for the metallic spectra was 
included in our determination. The fitting procedure is de- 
scribed in Sect.|4l This method relies on the ionisation and 
excitation balance in the star assuming LTE, but a number 
of potential problems can occur. Lines with any emission 
infilling must be avoided. Poor quality atomic data, undiag- 
nosed line blends and non-LTE effects may also be problem- 
atic. This also requires lines with a wide range of excitation 
energies and multiple ionisation states. Consequently, this 
technique could not successfully be applied to all stars in our 
sample. In a number of cases the excitation potential could 
not provide a useful constraint, but the ionisation potential 
was still useful. In these cases the constraint provided by 
Balmer line fits, combined with the ionisation balance could 
still provide fairly precise T^s and log g values. 

In order to place the stars of our sample on the H- 
R diagram, we needed to determine accurate luminosities. 
Photometric Johnson V magnitudes and B — V colours 
were derived fro m th e Hipp arcos V magnitudes using the 
conversion from IESAI (|l997l '). For stars without Hipparcos 
obse rvations, Joh nson magnitudes and colours were taken 
fromlVieira et a.1.1 (,200 3). Hcrbst & Shcvchcnko (1991), and 
Ide Winter et al. (|200l] ). Thes e magnitudes and colours are 
consistent with those used bv lAlecian et al.l (|2012al '). 

Distance were derived using Hipp a,rcos parallaxe s , whe n 
available, from the new reduction by Ivan LeeuwenI (|2007h . 
When no parallax was available, the star was associated 
with an OB association or star forming region, and a lit- 



erature dista nce to t hat wa s taken. Associations wer e taken 
from l Alecian et al.l (|2012al '). which agree with Vieira et al.l 
(|2003i r Uncertainties for these distances were taken to be 
the spatial extent of the associations, which were estimated 
from their angular size on the sky, assuming the associations 
are approximately spherical. For HD 169142 no association 
could reliabl y be made, and thus a literature photometric 
distance from lSvlvester et al.l (|l996l ) was used, and generous 
uncertainties were assumed. Literature sources for individ- 
ual distances, along with the distances themselves, can be 
found in Table, (3] 

Using these distances, bolometric luminosities were cal- 
c ulated. E(B — V) was calcula ted using the intrinsic colours 
of lKenvon fc HartmannI (|l995l ). and total extinction was cal- 
cul ated using an R y of 5 , as suggested for HAeBe stars 
bv [Hernandez et al.l (|2004 ). The bolometric correction of 
iBalonal ( 1994 ) was used. 

Using our luminosities and temperatures, we placed the 
stars on the H-R diagram, as shown in Fig. [T] Pre-main se- 
quence evolution ary tracks a nd isochrones were calculated 
with the CESAM (|More3 [l993) evolutionary code, assuming 
solar metallicities. The birth line, the locus of points on the 
H-R diagr am where a star becomes optically visible, was 
taken from iPalla fc StahleJ l| 19931 ). for an accretion rate of 
10~^ Mq yr~^. By comparison to these evolutionary tracks 
and isochrones, we determined stellar masses, and ages with 
respect to the birthline. These values, along with fractional 
pre-main sequence ages (r) and radii, are presented in Ta- 
ble [3l The uncertainties on mass and age for a star were 
based on the range of evolutionary tracks and isochrones 
that intersect the ellipse on the H-R diagram described by 
the uncertainties in the star's luminosity and Toff. Note that 
the choice of a birthline may introduce a further systematic 
u ncertainty into our ages. F or example, using the birthline 
of iBehrend fc Maedeij (|200ll ). with a mass dependent accre- 
tion rate, would generally increase our ages by ~0.5 Myr 
(but with the actual increase varying with mass). While we 
do not present uncertainties on r, the values should be con- 
sidered approximate. 



4 ABUNDANCE ANALYSIS 

The abundance analysis and determination of v sin i and 
microturbulence was performed by directly fitting synthetic 
spectra to the observations. Model spectra wer e calculated 
using the Zeeman sp ectrum synthesis code (|Landstreetl 
1 19881 : 1 Wade' et al.|[200ll ). This code performs polarised radia- 
tive transfer under the assumption of LTE. Optimisations to 
the code for stars with negligible magnetic fields have been 
included. A Levenber g-Marquardt minimisation proce- 
dure (e.g. iPress et al.lil992l ) was used to determine best fit 
parameters. 

The optimisations to the Zeeman spectrum synthesis 
routine exploit symmetries in a non-magnetic chemically 
homogeneous model of a star to dramatically reduce the 
amount of computation required. As in the original version 
of Zeeman, the visible disk of the star is divided into a 
number of surface elements. In this version it is assumed 
that the local emergent spectra only vary with projected ra- 
dial distance from the centre of the disk of the star. Thus 
radiative transfer is only performed for a set of surface el- 
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Figure 1. H-R diagram for tlie stars in tliis study. Evolutionary tracks (solid lines) are labelled by mass in solar masses. Isochrones 
(dashed lines) are labelled by age in Myr. The birth line (right dot-dashed line) for an accretion rate of 10~^ Mq yr~^, and the zero-age 
main sequence line (ZAMS, left dot-dashed line) are also shown. Circles are chemically normal stars, squares are A Boo stars, the diamond 
is the possible Bp star V380 Ori A. Open symbols are stars with confirmed magnetic field detections. 



ements with different radial distances, and the results are 
reused for different angular positions, which minimises the 
number of times radiative transfer must be performed. The 
optimisations also assume the local line absorption (Voigt) 
and anomalous dispersion (Faraday- Voigt) profiles only vary 
vertically through the atmosphere, and not with position on 
the stellar disk, substantially reducing the number of times 
these profiles must be calculated. In this version of the code, 
limb darkening is still calculated directly, by performing ra- 
diative transfer at different radial positions on the disk of 
the star. Rotational broadening is also calculated directly 
by summing local emergent spectra that have been Doppler 
shifted across the disk of the star. 

Detailed comparisons of the optimised code to the origi- 
nal version have been performed. For the same input param- 
eters, with no magnetic field, identical spectra are produced. 
The original code has been ch ecked against oth er spec- 
trum synthesis codes in detail bv lWade et al.l (|200ll ). with a 
very good agreement found. Similarly, the optim ised version 
of Ze eman has been checked against SynthS (jKochukhovl 
[2003) and again a good agreement was found. 

The Levenber g-Marquardt minimisation procedure 
JPress et al ] |1992') provides a robust and efficient routine 
for iteratively fitting synthetic spectra to observations. Best 
fit models from this routine have been compared to fits ob- 



tained by hand for a wide variety of stars, and good agree- 
ment between the two methods has been consistently found. 

Input atomic data w ere taken from the Vienna Atomic 
Line Database (VALD) (|Kupka et al.lll999l '). using an 'ex- 
tract stellar' request and solar chemical abundances. Up- 
dated line lists with peculiar abundances were extracted as 
necessary. Model atmosph eres were computed using the AT- 
LAS9 code l|Kurucdll99^ ) with solar abundances. Model at- 
mospheres were calculated on a grid in steps of 250 K in Tcb 
and 0.5 in log g, then interpolated linearly when more pre- 
cise values were desired. Interpolation on this grid of models 
introduces less then 1% relative error to the line depths in 
the resulting spectrum. 

Spectra were initially fit simultaneously for chemical 
abundances, i;sini, microturbulence, and radial velocity, us- 
ing the best fit T^g and log g determined from the Balmer 
lines. Then, if a good constraint could be obtained, the fits 
were repeated including Tcb and log g as free parameters. If 
this produced a good fit with well constrained parameters, 
then these TcS and log g values were used instead of the 
initial Balmer line fit values. If only one of T^s and log g 
was well constrained from the metallic line fit, then the fit 
was repeated using the Balmer line value for the uncon- 
strained parameter. We considered a parameter to be well 
constrained if the fitting routine reliably converged to the 
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Table 3. Fundamental parameters derived for the stars in this study, t is the fractional pre-main sequence age. The 'Magnetic' column 
lists which star s have confirmed m agnetic field detections (M), and w h ich stars we consid er to be non-magnetic (N ) . The references for 
distances are: "Ivan Leeuwenll2007l (Hipparcos); Ide Zeeuw et al.lll999l: terown et al Ill994 toolan fc MathieullioOll : HVieira et alll2003l: 
■ ISvlvester et al.lll996l . The luminositv for HD 190073 was taken from lCatala et al.l \20ff^ (indicated by *). 



ID 


Toff (K) 




log g 




distance (pc) 


log(L/L0) 


M/Mq 


R/Rq 


age (Myr) 


r 


Magnetic 


HD 17081 


12900 ± 400 


3, 


,8±0, 


,2 


120 ±3" 


2.67 ±0.04 


4.4 ±0.2 


4.3 ±0.3 


'^•'^-0.1 


0.48 


N 


HD 31293 


9800 ± 700 


3, 


,9±0, 


,3 


139 ± 19'' 


1.77 ±0.13 


2.55 ±0.2 


2.7±0.5 


^■"-0.6 


0.65 


N 


HD 31648 


8800 ± 190 


4, 


.1 ±0, 


.2 


137 ± 25'' 


1.27 ±0.16 


2.1 ±0.25 


1.9 ±0.4 


"■^-2.0 


0.84 


N 


HD 36112 


8190 ±150 


4, 


.1 ±0, 


.4 


279 ± 70" 


1.82 ±0.22 


2.8 ±0.5 


4.0 ±1.0 


^^■''-1.0 


0.52 


N 


HD 68695 


9000 ± 300 


4, 


,3±0, 


,3 


410 ±36'' 


1.34 ±0.08 


2.2 ±0.15 


1.9 ±0.2 


c 9+3.0 

°-^-i.o 


0.76 


N 


HD 139614 


7600 ± 300 


3, 


,9±0, 


,3 


147 ± 37'' 


0.88 ±0.22 


1.7±0.1 


1.6 ±0.4 




0.82 


N 


HD 141569 


9800 ± 500 


4, 


,2±0, 


,4 


116 ± 8° 


1.52 ±0.07 


2.4 ±0.2 


2.0 ±0.3 


0.83 


N 


HD 142666 


7500 ± 200 


3, 


.9±0, 


.3 


145 ± 18'' 


1.28 ±0.11 


1.95 ±0.15 


2.5 ±0.3 




0.63 


N 


HD 144432 


7400 ± 200 


3, 


.9±0, 


.3 


160 ± 29'^ 


1.27 ±0.16 


1.95 ±0.2 


2.6 ±0.5 


5 5+2-5 


0.58 


N 


HD 163296 


9200 ± 300 


4, 


,2±0, 


,3 


119 ± 11" 


1.49 ±0.08 


2.3 ±0.1 


2.2 ±0.2 


4 1+1-0 


0.71 


N 


HD 169142 


7500 ± 200 


4, 


,3±0, 


,2 


145 ± 50-f 


1.01 ±0.30 


1.7±0.2 


1.9 ±0.7 


^■^-3.0 


0.52 


N 


HD 176386 


11000 ± 400 


4, 


,1 ±0, 


,3 


128 ± 13° 


1.79 ±0.09 


2.8 ±0.2 


2.2 ±0.3 


9 7+1.0 
^- '-1.0 


0.93 


N 


HD 179218 


9640 ± 250 


3, 


.9±0, 


.2 


254 ± 38" 


2.04 ±0.13 


3.1 ±0.3 


3.7 ±0.6 


1 1+0-'^ 


0.55 


N 


HD 244604 


8700 ± 220 


4, 


.0±0, 


.2 


380 ± 79= 


1.83 ±0.18 


2.75 ±0.4 


3.6 ±0.8 


1 9+^-1 
1^-^-1.0 


0.61 


N 


HD 245185 


9500 ± 750 


4, 


.0±0, 


.4 


450 ± 50'* 


1.44 ±0.11 


2.3 ±0.2 


1.9 ±0.4 


5 5+2-" 


0.90 


N 


HD 278937 


8000 ± 250 


4, 


,1 ±0, 


,2 


318 ±43'' 


1.04 ±0.12 


1.8 ±0.1 


1.7±0.3 


9 5+5-° 


0.76 


N 


T Ori 


8500 ± 300 


4, 


,2±0, 


,3 


380 ± 44'= 


1.67 ±0.10 


2.45 ±0.15 


3.1 ±0.4 


Q Q + 1.2 

'-'■'-'-1.3 


0.75 


N 


HD 101412 


8600 ± 300 


4, 


,0±0, 


,5 


600 ± 100*= 


1.92 ±0.14 


3.0 ±0.3 


4.2 ±0.8 


1 9+0.8 
l--^-0.7 


0.50 


M 


HD 190073 


9230 ± 260 


3, 


.7±0, 


.3 


> 340° 


1.92 ±0.12* 


2.9 ±0.3 


3.6 ±0.5 


1 6+°-'' 

l^-"-0.6 


0.62 


M 


V380 Ori A 


12600 ± 1000 


4, 


.0±0, 


.5 


398 ±91'= 


2.66 ±0.20 


4.4 ±0.7 


4.5 ±1.2 


0.01-5 


0.55 


M 


V380 Ori B 


5800 ± 350 


4, 


.1 ±0, 


.3 


398 ±91'= 


1.62 ±0.27 


3.3 ±0.8 


6.4 ±2.2 


0.0 


N 



same value for diflterent initial conditions, and if that value 
was consistent with the Balmer line profiles. The fits were 
generally performed on 5 independent regions ~500 A long. 
The approximate wavelength ranges usually were 4400-4800, 
4900-5500, 5500-6000, 6000-6500, and 6600-7600 A, and oc- 
casionally 4150-4280 A, with significant gaps due to non- 
photospheric features. These windows typically contained 
several hundred spectral lines (the number varying greatly 
with wavelength and Tofj), however many of these lines are 
very weak and hence do not provide significant constraints 
on the best fit abundances or atmospheric parameters. The 
exact wavelength ranges varied between stars, due to vary- 
ing regions of emission and different Balmer line widths. 

Microturbulence and v sin i were determined by min- 
imisation, simultaneously with the other stellar parameters, 
using the entire spectral ~500 A window. This method relies 
on resolved, rotationally broadened spectral lines for con- 
straining vsini, and a range of both weak and strong lines 
for constraining microturbulence. 

We verified this method of determining Tcff and log g 
by analysing well studied stars, and comparing to literature 
results. In this study we have analysed tt Get, and compared 
our re s ults with the very precise study by iFossati et al.l 
(|2009l '). iFossati et all \200^ ) used photometry and spectral 
energy distributions, as well as Balmer lines, and excitation 
and ionisation balances when determining their parameters. 
They also used different software tools to perform their mod- 
elling, thus their results are truly independent. Our results 
are fully consistent with theirs, as discussed in Sect. lAll 
which demonstrates the accuracy of our methodology. 

For this method to be successful, lines contaminated by 



emission must be avoided. While a line entirely in emission is 
easily identified, lines with small amounts of emission infill- 
ing are harder to identify. When multiple observations were 
available, variability in the emission could be used to iden- 
tify lines with emission infilling. If only one observation was 
available, special attention was paid to lines with low ex- 
citation potentials and large oscillator strengths, such lines 
being more likely to contain emission. If inconsistent fits 
were obtained between normal and low excitation potential 
lines, and no other clear explanation for the inconsistency 
could be found (such as an error in Tcff or the atomic data), 
the low excitation potential lines were considered to likely 
contain emission and were excluded from the final fit. Ad- 
ditionally, the shapes of line profiles were examined and, if 
they could not be reproduced by the synthetic spectra, they 
were considered to be contaminated by a non-photospheric 
component. The presence of veiling in our spectra is assessed 
in detail in Section 16.41 and we find it is not present at a 
significant level. 

The use of multiple spectral windows is valuable, as 
it allows us to take the average and standard deviation of 
abundances determined across the windows. The standard 
deviation in particular provides a reasonable estimate of the 
uncertainty of derived parameters. Since line-to-line abun- 
dances scatter is introduced by errors in the atomic data 
as well as errors in the model atmosphere, these effects are 
included in the standard deviation. The standard deviation 
may underestimate the impact of such errors, but provides a 
robust and easily understandable uncertainty estimate. The 
use of multiple spectral windows also allows us to verify 
the atmospheric parameters derived from metal line fits. If 
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the parameters difTer substantially from window to window, 
they are likely poorly constrained in some if not all windows, 
in which case we can rely more heavily on the parameters 
from Balmer line fitting. 

For chemical elements with fewer than four lines provid- 
ing good constraints on the abundance in any spectral win- 
dow (including all cases of elements with lines in only one or 
two spectral windows) the uncertainties were estimated by 
eye rather than using a standard deviation. These uncertain- 
ties were chosen to include line-to- line scatter, uncertainties 
in the atmospheric parameters, and potential normalisation 
errors. 

Examples of typical best fit synthetic spectra compared 
to the observed spectrum of HD 139614 in the 5000-5050 A 
and 6100-6200 A windows are presented in Fig. [21 Final av- 
eraged best fit abundances and atmospheric parameters are 
presented in Table |4l together with uncertainties. Elements 
with abundances based on less than four lines are indicated 
with an asterisk in the table. 

We find 11 stars with varying degrees of A Boo peculiar- 
ity, one star (V380 Ori A) that appears to be a weak Ap/Bp 
star, and 9 stars that are chemically normal (although re- 
sults for V380 Ori B are very uncertain). The properties of 
the individual stars are discussed in detail in the Appendix. 
The final best fit abundances are plot ted relative to solar 
abundances from iGrevesse et all (|2005l ) in Fi gures \AZl IA31 
IA4I and lASI Illustrations of our fits for all of the stars in our 
study can be found in Figures [A6I to [A25l 



5 TRENDS 

In order to investigate trends in chemical abundance rela- 
tive to other stellar parameters, we constructed a peculiarity 
index to describe whether a star is chemically normal, more 
like a A Boo star, or more like an Ap/Bp star. This index 
is designed to be positive for stars with overabundances of 
iron peak elements, zero for chemically normal stars, and 
negative for stars underabundant in iron peak elements. The 
peculiarity index is the difference between the average of the 
Cr, Fe, and Ni abundances relative to solar and the average 
of the C, N, and O abundances relative to solar. Specifically, 
the index is calculated as: 

[P] = i(([Cr] + [Fe] + [Ni]) - ([C] + [N] + [O])), (1) 

where the abundances are relative to solar: [X] = 
\og{Nx,*/Ntot,*) - \og{Nx,e/Ntot,e)- These particular ele- 
ments were chosen because they tend to have peculiar abun- 
dances in chemically peculiar stars (for Cr, Fe, and Ni) or 
normal abundances in chemically normal stars (for C, N and 
O). These elements were also chosen because reliable abun- 
dances for them are available for most stars in our sample. In 
the few cases in which we could not determine an abundance 
for an element used in the peculiarity index, that element 
was left out of the average. 

We plotted the index [P] against effective temperature, 
mass, age, and fractional pre-main sequence age, as shown 
in Fig. [3] No clear trends were found, with a large scatter 
in chemical abundance at most temperatures, masses, and 
ages. This peculiarity index was also compared to stellar ra- 
dius, Dsini, and microturbulence, but again no clear trends 
were found. Thus we conclude that the A Boo peculiarities 



that we see are not restricted to a narrow mass range, and 
that they are present throughout the pre-main sequence. 

6 DISCUSSION 

6.1 A Boo peculiarities 

The most striking result of this study is the discovery of 
a large number of HAeBe stars with moderate A Boo type 
chemical peculiarities. This trend seems to be present in 
both stars with and without detected magnetic fields, as 
shown by HD 101412. We find 7 stars showing fairly strong 
A Boo peculiarities, and another 4 stars with weaker A Boo 
peculiarities, while 8 stars are chemically normal and 1 is 
rather uncertain but appears to be normal. With the de- 
rived abundances, we see a continuous distribution of A Boo 
peculiarities, from fairly strong to nearly undetectable. This 
incidence of roughly 50% is much higher than the incidence 
of A B oo stars on the main sequence of roug hly 2% l|Paunzenl 
l200lh . While a detection rate of 50% may be an overestimate 
(since our sample is biased towards stars with limited emis- 
sion in their optical spectra), the true incidence of A Boo 
peculiarities in HAeBe stars certainly is much higher than 
observed on the main sequence, suggesting an important 
connection with the pre-main sequence evolutionary phase. 

The sample of stars in our study is biased towards stars 
with weaker circumstellar emission or absorption in their 
spectra, since these stars allow for a more complete anal- 
ysis of their photospheric spec tra. However, t he star s were 
drawn from the larger sample of lAlecian et all l|2012al ) which 
did not have this bias. We found 65% of the stars in their 
sample had modest amounts of circumstellar contamination, 
and would have been suitable for our photospheric analysis. 
Thus, as a lower limit, 33% of HAeBe stars should show A 
Boo peculiari ties, assuming that al l the unsuitable stars in 
the sample of lAlecian et al.l ((2012^ are chemically normal, 
and using our derived incidence of peculiarities for the suit- 
able stars. The true incidence is likely higher than this lower 
limit. 

There are a lmost certainly more u ndetected A Boo stars 
in the sample of lAlecian et al.l (|2012al ) that we did not anal- 
yse, simply because analysing their com plete sample was be - 
yond the scope of this study. Indeed lAlecian et al.l (|2012ar ) 
mention that HD 34282 appears to have weak lines of iron 
peak elements, but a normal O line. To roughly estimate the 
number of potential A Boo stars in the remaining sample, 
we compared synthetic spectra computed using literature 
values for their atmospheric parameters and assuming solar 
abundances with the ESPaDOnS observations. In doing so 
we identified 5 additional stars that have suspiciously weak 
metal lines for their literature temperatures: HD 37806, HD 
98922, HD 144668, HD 76534, and VX Cas. If all these stars 
do indeed show A Boo pe culiarities, 35% of the stars ob- 
served with ESPaDOnS bv lAlecian et all (12012^ 1 would be 
A Boo stars, which is consistent with our results. 

This large incidence of A Boo type peculiarities in 
HAeBe stars is qualitatively consistent with a selective ac- 
cretion model of A Boo stars. HAeBe stars have recently 
undergone accretion, and may still be accreting, so it is rea- 
sonable that a process which depends on accretion is seen 
more frequently in these stars. Thus we consider our results 
supportive of the selective accretion model of A Boo stars. 
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Table 4. Best fit parameters for the stars in our sample. Cliemical abundances are in units of log{N x / Ntot) ■ Elements marked with an 
asterisk are based on less then ~4 useful lines and have uncertainties estimated by eye. Microturbulence is given by ^. 





HD 17081 


HD 31293 


HD 31648 


HD 36112 


HD 68695 


HD 139614 


HD 141569 


Solar 




(tt Oct) 


(AB Aur) 














TLtf (K) 


12900 ± 400 


9800 ± 700 


8800 it 190 


8190 ± 150 


9000 it 300 


7600 it 300 


9800 it 500 




log g (cgs) 


3.8 ± 0.2 


3.9 ± 0.3 


4.1 it 0.2 


4.1 it 0.4 


4.3 it 0.3 


3.9 ± 0.3 


4.2 it 0.4 




V sin i (km ) 


20.9 ± 1.2 


116 ± 9 


101.2 it 1.7 


57.8 it 1.0 


51 it 4 


25.6 it 0.4 


222 it 7 




^ (kms~l) 


1.7 ± 1.0 


^ 4 


3.2 ± 1.1 


2.97 it 0.24 


1.3 ± 1.1 


3.7 it 0.4 


^ 2 




He 


-0.91 ± 0.14 


-1.24 ± 0.20* 










-1.21 it 0.37* 


-1.11 


C 


^ -3.4* 


-3.33 ± 0.22 


-3.67 it 0.07 


-3.61 it 0.16 


-3.27 it 0.18 


-3.75 ± 0.22 


-3.63 it 0.29 


-3.65 


N 


-4.03 ± 0.15* 


-3.9 ± 0.3* 






-3.7 ± 0.3* 


-4.17 ± 0.13 


-4.0 it 0.3* 


-4.26 


O 


-3.16 ± 0.13 


-3.27 ± 0.20 


-3.28 ± 0.05* 


-3.18 ± 0.10* 


-3.17 it 0.10 


-3.33 it 0.10* 


-3.05 it 0.10 


-3.38 


Ne 


-3.76 ± 0.28 














-4.20 


Na 


-5.30 ± 0.10* 




-5.72 ± 0.25* 


-5.58 ± 0.15* 


-6.2 ± 0.4* 


-6.14 ± 0.12 


^ -5.2* 


-5.87 


Mg 


-4.43 ± 0.15 


-4.89 ± 0.33 


-4.12 it 0.13 


-4.23 it 0.14 


-5.11 it 0.14 


-4.74 ± 0.07 


-4.90 ± 0.11 


-4.51 


Al 


-5.81 ± 0.17 












-5.6 ± 0.5* 


-5.67 


Si 


-4.55 ± 0.14 


-4.76 ± 0.22 


-4.29 it 0.19 


-4.58 it 0.11 


-5.36 it 0.15 


-5.10 ± 0.13 


-4.99 it 0.31 


-4.53 


P 


-6.43 ± 0.14 














-6.68 


S 


-4.99 ± 0.08 




-4.6 ± 0.3* 


-4.79 ± 0.16* 


-4.6 ± 0.3* 


-5.18 ± 0.15 




-4.90 


Ar 


-5.6 ± 0.3* 














-5.86 


K 












-7.29 it 0.25* 




-6.96 


Ca 


-5.7 ± 0.2* 


-6.10 ± 0.17 


-5.38 ± 0.23 


-5.59 ± 0.16 


-6.47 it 0.23 


-6.18 ± 0.16 


-6.34 it 0.38 


-5.73 


Sc 




-9.04 ± 0.34 


-8.72 ± 0.15 


-8.90 ± 0.26 


-9.49 ± 0.26 


-9.39 ± 0.15 


-9.18 ± 0.37 


-8.87 


Ti 


-7.37 ± 0.15 


-7.47 ± 0.27 


-6.82 it 0.09 


-6.98 ± 0.20 


-7.71 it 0.25 


-7.51 ± 0.14 


-7.74 it 0.32 


-7.14 


V 








-7.8 ± 0.4* 




-8.4 it 0.5* 




-8.04 


Cr 


— D.oi in u.J-O 




fi nQ 4-017 


— D. ou HI u. r 1 


QS 4- n OA 
— o.yo HI yj.^'i 


— O.OU HI U.J- J- 


7 '^1 4- n ^1 

— i .OJ- HI u.or 


Q 


Mn 


-6.6 ± 0.3* 




-6.59 it 0.08 


-6.70 it 0.16 


-6.4 ± 0.5* 


-7.43 ± 0.19 




-6.65 




— 4. ± U.OK 


— 4.yi it O.^O 


— 4.4i itO.lO 


A ACi -I- r\ 1 A 
— 4.4y zt U.i4 


— 0. iO It u.zz 


K A7 r\ -t ^ 

— o.O / ± U. lo 


OC -LA 
— O.ZO It O.OZ 


— 4.59 


Co 
















-7.12 


Ni 


-5.85 ± 0.06 


-6.16 ± 0.27 


-5.66 ± 0.05 


-5.76 it 0.20 


-6.25 ± 0.32 


-6.33 ± 0.14 


sS -6.0* 


-5.81 


Cu 








!j -7.5* 




-8.6 it 0.4* 




-7.83 


Zn 








-7.8 ± 0.4* 




-8.3 it 0.3* 




-7.44 


Sr 




—9.6 ± 0.3* 




— 9.4 ± 0.3* 






-10.5 it 0.5* 


—9.12 


Y 






-9.57 it 0.30 


-9.56 ± 0.17 




-10.23 it 0.15 




-9.83 


Zr 








-9.3 ± 0.3* 








—9.48 


Ba 




—9.7 ± 0.4* 


-9.56 ± 0.24* 


—9.46 it 0.29 


— 10.3 it 0.4* 


-10.29 it 0.21 


-9.2 it 0.8* 


—9.87 


















-10.91 


Co 












^ —10.5* 




-10.34 


Nd 
















— 10.59 


Eu 
















— 11.52 




HD 142666 


HD 144432 


HD 163296 


HD 169142 


HD 176386 


HD 179218 


HD 244604 


Solar 


TLff (K) 


7500 ± 200 


7400 it 200 


9200 it 300 


7500 it 200 


11000 it 400 


9640 ± 250 


8700 ± 220 




log g (cgs) 


3.9 ± 0.3 


3.9 ± 0.3 


4.2 it 0.3 


4.3 it 0.2 


4.1 it 0.3 


3.9 ± 0.2 


4.0 it 0.2 




V sin i (km s^ ) 


68.0 ± 0.9 


80.3 ± 1.0 


122 ± 3 


51.6 ± 0.5 


169.0 ± 1.5 


70 it 4 


101 it 5 




5 (kms^^) 


3.55 ± 0.31 


3.62 ± 0.23 


1.5 it 1.2 


2.09 it 0.47 


1.7 it 0.7 


2.0 ± 0.5 


1.9 it 0.4 




He 










-1.17 ± 0.20* 


-1.15 it 0.15* 




-1.11 


C 


-3.62 ± 0.15 


-3.79 ± 0.18 


-3.82 it 0.25 


-3.55 it 0.12 


-3.38 ± 0.29* 


-3.45 ± 0.16 


-3.69 it 0.17 


-3.65 


N 










-3.9 ± 0.3* 


-3.5 it 0.4* 


-4.2 it 0.3* 


-4.26 


O 


-3.18 ± 0.15* 


-3.17 ± 0.10* 


-3.31 it 0.15 


-3.38 it 0.13 


-3.20 ± 0.10* 


-3.10 ± 0.13 


-3.23 it 0.08 


-3.38 


Ne 
















-4.20 


Na 


-5.77 ± 0.15* 


-5.86 ± 0.09 


-5.6 ± 0.5* 


-6.18 ± 0.08* 




-5.5 it 0.3* 


-5.4 it 0.3* 


-5.87 


Mg 


-4.61 ± 0.06 


-4.46 it 0.04 


-4.12 ± 0.15 


-4.89 it 0.06 


-4.28 ± 0.15 


-4.73 it 0.17 


-4.03 ± 0.22 


-4.51 


Al 






-5.5 it 0.4* 




-5.5 it 0.3* 


-6.0 ± 0.4* 




-5.67 


Si 


-4.91 ± 0.20 


-4.81 it 0.24 


-4.24 ± 0.14 


-5.13 ± 0.10 


-4.44 ± 0.14 


-4.82 it 0.22 


-4.34 ± 0.06 


-4.53 


P 
S 


-4.70 ± 0.15 


-4.82 it 0.05 




-5.10 ± 0.12 




-4.2 ± 0.4* 


-4.48 ± 0.15* 


-6.68 
-4.90 


Ar 
















-5.86 


K 
















-6.96 


Ca 


-5.98 ± 0.18 


-5.78 ± 0.15 


-5.54 it 0.30 


-6.15 it 0.12 


-6.27 ± 0.19* 


-6.25 ± 0.38 


-5.42 it 0.24 


-5.73 


Sc 


-9.23 ± 0.23 


-9.09 it 0.15 


-8.76 ± 0.08 


-9.58 ± 0.24 


-9.0 ±0.4* 


-9.40 ± 0.16 


-8.75 ± 0.24 


-8.87 


Ti 


-7.39 ± 0.20 


-7.26 it 0.16 


-6.92 it 0.10 


-7.65 ± 0.10 


-7.12 ± 0.29 


-7.57 it 0.12 


-6.81 it 0.31 


-7.14 


V 


-8.4 ± 0.4* 


-8.3 ± 0.4* 




-8.49 ± 0.16* 








-8.04 


Cr 


— D.DU ± 


— D.Ol It 0.^3 


K AK -X- A OA 

— o.yo ± O.ZO 


no -La 1 q 
— b.yo It 0. lo 


— D.4D ± O.Zl 


f; A7 -1- A 11 
— D.y / It U. iO 


— D.O/ itO.lo 


— 6.40 


Mn 


-7.08 ± 0.15 


-6.96 it 0.15 


-7.0 it 0.5* 


-7.57 it 0.20 






-6.91 it 0.23 


-6.65 


Pc 


-4.84 ± 0.11 


-4.70 it 0.09 


-4.39 ± 0.15 


-5.13 it 0.11 


-4.48 ± 0.27 


-5.03 it 0.13 


-4.35 ± 0.24 


-4.59 


Co 


-6.8 ± 0.4* 


-7.1 ± 0.5* 




^ -7.6* 








-7.12 


Ni 


-6.16 ± 0.13 


-6.00 ± 0.15 


-5.72 ± 0.12 


-6.40 ± 0.16 




^ -5.9* 


-6.00 it 0.18 


-5.81 


Cu 




< -7.7* 




-8.85 ± 0.28 








-7.83 


Zn 


-7.9 ± 0.3* 


-7.57 it 0.20 




-8.71 ± 0.08* 








-7.44 


Sr 


-9.56 ± 0.19* 


-9.4 it 0.3* 




-9.58 ± 0.11* 




-10.2 ± 0.2* 




-9.12 


Y 


-9.94 ± 0.22 


-9.80 ± 0.13 




-10.18 ± 0.14 




-10.0 ± 0.4* 


-9.8 it 0.4* 


-9.83 


Zr 


-9.36 ± 0.26* 


-9.3 it 0.4* 




-9.65 ± 0.13* 








-9.48 


Ba 


-9.81 ± 0.18 


-9.69 it 0.16 


-10.0 it 0.4* 


-9.96 it 0.21 


< -9.0* 


-10.4 it 0.4* 


-9.6 it 0.2* 


-9.87 


La 


^ -10.8* 


-11.2 it 0.4* 












-10.91 


Ce 
















-10.34 


Nd 








-10.64 ± 0.14* 








-10.59 


Eu 








-12.19 it 0.10* 








-11.52 
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Table 4 - continued Best fit parameters for the stars in our sample. Chemical abundances are in units of log{N x / Ntot) ■ Elements 
marked with an asterisk are based on less then ~4 useful lines and have uncertainties estimated by eye. Microturbulence is given by ^. 





HD 245185 


HD 278937 
(IP Per) 


T Ori 
(BD -05 1329) 


HD 101412 


HD 190073 


V380 Ori A 


V380 Ori B 


Solar 


Td! (K) 


9500 ± 750 


8000 ± 250 


8500 ± 300 


8600 ± 300 


9230 ± 260 


12600 ± 1000 


5800 ± 350 




log g (cgs) 


4.0 ± 0.4 


4.1 ± 0.2 


4.2 ± 0.3 


4.0 ± 0.5 


3.7 ± 0.3 


4.0 ± 0.5 


4.1 ± 0.3 




V sin i (km ) 


136 ± 10 


83.8 ± 4.6 


163 ± 11 


6.8 ± 0.4 


8.50 ± 0.23 


9.9 ± 1.0 


24.7 ± 1.9 




5 (kms~^) 


^ 4 


2.0 ± 0.9 


2.6 ±1.1 


^ 2 


^ 2 


^ 3 


^ 2 




Ho 


-0.9 ±0.3* 




-1.4 ± 0.5* 


-0.9 ± 0.3* 


-0.90 ± 0.20* 


-0.67 ± 0.20* 




-1.11 


C 


-3.72 ± 0.15 


-3.57 ± 0.17 


-3.59 ± 0.19 


-3.58 ± 0.12 


-3.72 ± 0.14 


-3.4 ± 0.4* 




-3.65 


N 


-3.8 ± 0.4* 


-4.0 ± 0.3* 


-4.2 ± 0.4* 




-3.60 ± 0.20* 






-4.26 


O 


-3.17 ± 0.17 


-3.41 ± 0.05 


-3.16 ± 0.12 


-3.12 ± 0.09 


-3.26 ± 0.11 


-3.00 ± 0.16 




-3.38 


Nc 










-3.8 ± 0.3* 


-3.93 ± 0.36 




-4.20 


Na 




-6.3 ± 0.3* 


-6.1 ± 0.3* 


-5.56 ± 0.15* 


-5.13 ± 0.15* 




-5.95 ± 0.30* 


-5.87 


Mg 


-5.39 ± 0.23 


-4.91 ± 0.11 


-5.05 ± 0.23 


-4.95 ± 0.15 


-4.38 ± 0.17 


-4.19 ± 0.14 


-4.7 ± 0.7* 


-4.51 


Al 


-5.4 ± 0.5* 




-5.2 ± 0.4* 


-5.99 ± 0.20* 


-5.58 ± 0.17* 


-5.07 ± 0.21* 




-5.67 


Si 


-4.95 ± 0.16 


-4.97 ± 0.19 


-4.89 ± 0.20 


-5.30 ± 0.19 


-4.48 ± 0.06 


-4.11 ± 0.09 


-4.44 ± 0.26 


-4.53 


P 












-5.86 ± 0.30* 




-6.68 


S 




-4.83 ± 0.14 


-4.1 ± 0.3* 


-4.96 ± 0.12 


-4.57 ± 0.20* 


-4.93 ± 0.26 




-4.90 


Ar 
















-5.86 


K 
















-6.96 


Ca 


-6.12 ± 0.25 


-6.28 ± 0.18 


-6.00 ± 0.31 


-6.37 ± 0.20 


-5.61 ± 0.12 


s; -5.6* 


-5.78 ± 0.40 


-5.73 


Sc 


-9.6 ± 0.7* 


-9.52 ± 0.21 


-9.4 ± 0.3* 


-9.33 ± 0.17 


-8.85 ± 0.20 






-8.87 


Ti 


-7.93 ± 0.38 


-7.81 ± 0.10 


-7.55 ± 0.17 


-7.75 ± 0.23 


-7.16 ± 0.12 




-6.78 ± 0.28 


-7.14 


V 


^ -7.5* 






-8.1 ± 0.3* 


-7.99 ± 0.25* 




-7.50 ± 0.79 


-8.04 


Cr 


-7.27 ± 0.34 


-7.04 ± 0.18 


-6.82 ± 0.07 


-6.69 ± 0.29 


-6.13 ± 0.13 


^ -4.8* 


-5.92 ± 0.26 


-6.40 


Mn 


-7.0 ± 0.6* 


-7.3 ± 0.4* 


^ -6.0* 


-6.98 ± 0.09 


-6.54 ± 0.08 


-5.92 ± 0.16 


-5.97 ± 0.43 


-6.65 


Fc 


-5.27 ± 0.33 


-5.12 ± 0.16 


-4.98 ± 0.10 


-5.08 ± 0.19 


-4.42 ± 0.06 


-3.92 ± 0.12 


-4.25 ± 0.16 


-4.59 


Co 














-6.6 ± 0.6* 


-7.12 


Ni 


-6.5 ± 0.7* 


-6.37 ± 0.22 




-6.13 ± 0.27 


-5.69 ± 0.17 


-5.19 ± 0.18 


-5.58 ± 0.32 


-5.81 


Cu 




^ -7.5* 




-8.6 ± 0.3* 






-7.4 ± 0.7* 


-7.83 


Zn 




^ -8.0* 






-7.52 ± 0.25* 




-7.4 ± 0.6* 


-7.44 


Sr 






-9.3 ± 0.2* 










-9.12 


Y 




-10.4 ± 0.4* 




-10.12 ± 0.25 


-9.70 ± 0.20 




^ -8.2* 


-9.83 


Zr 








-9.8 ± 0.3* 


-9.2 ± 0.3* 






-9.48 


Ba 


-10.0 ± 0.7* 


-10.2 ± 0.3* 


-10.3 ± 0.3* 


-10.71 ± 0.24 


-9.71 ± 0.13* 




-8.89 ± 0.40 


-9.87 


La 
















-10.91 


Co 










^ -9.6* 




^ -8.5* 


-10.34 


Nd 










-10.2 ± 0.3* 






-10.59 



The large number of pre-main sequence A Boo stars we 
find suggests that the majority of A Boo stars may develop 
their peculiarities on the pre-main sequence by accreting 
their circumstellar material. However, we cannot rule out 
the existence more of evolved A Boo stars that only develop 
peculiarities on the main sequence, by accreting from dif- 
fuse interstellar clouds. Due to the relatively short lifetime 
of A Boo peculiarities after accretion has halted, the pres- 
ence of such peculia rities later in the main sequence (e.g. 
iPaunzen et al.l I2OO2I ') may require a different source of ac- 
creted material from that of pre-main sequence A Boo stars. 

iTurcotte fc Charbonneaul (|l993l ) modelled the diffusion 
of chemical elements in the atmosphere of a star (with Tcff — 
8000 K and log g = 4.3) that was accreting gas depleted in 
iron peak elements. They found that A Boo peculiarities can 
be generated quickly (~0.1 Myr), but also dissipate quickly 
after accretion has stopped (~1 Myr). This implies that, in 
the context of selective accretion, there is a good chance that 
HAeBe stars with A Boo peculiarities are still accreting, or 
were very recently accreting. 

The range in strengths of the A Boo peculiarities seen 
in the stars in our sample might be interpreted as the con- 
sequence of different accretion rates in the stars. However, 
iTurcotte fc Charbonneaul l| 19931 ) do not find a strong impact 
on surface abundances from different accretion rates, as long 
as they are sufficient to overwhelm diffusion and rotational 
mixing. 



Alternately, the differences in the strengths of peculiar- 
ities could reflect the abundances of the accreted gas. These 
differences being due, for example, to the efficiency with 
which iron peak elements are bound into dust grains, or to 
the degree to which dust is blown away from the star while 
gas is accreted. 

In this context, the chemically normal stars are most 
likely not accreting significant amounts of material, and thus 
any A Boo peculiarities have dissipated. However, it is pos- 
sible that accretion is ongoing but the selection process in 
the circumstellar material has broken down. For example, 
dust grain formation may not be proceeding efficiently, and 
thus the accreted gas would not be depleted in iron peak 
elements 

Unfortunately, deriving accretion rates for HAeBe stars 
is not a simple process. For example, Donohcw fc BrittainI 
(|201ll ) measured accretion rates based on Balmer jump 
excess fluxes. They obtained values between a few 10~^ 
and 10~^ M(r,yT~^. T hese values rely on models from 
iMuzeroUe et al.! (|2004h . which are constructed assuming 
magnetospheric accretion. Given the absence of detectable 
magnetic flelds in most of the stars in our study, mag- 
netospheric accre tion may not be an accurate model (see 
IWade et al.ll2007l ). and thus these accretion rates may not 
be accurate either. Generally, accretion rates derived for 
HAeBe stars are based on rather simple models, and thus 
it is not clear how reliable those accretion rates are. Eleven 
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Figure 2. Comparison of the observed spectrum (jagged line) to the best fit synthetic spectrum (smooth line) for HD 139614. Two 
independent wavelength regions are presented. Lines have been labelled by their major contributing species. 



of the stars in our sample h ave measured accretion rates 
bv lDonehew fc BrittainI l|201 lY Supplementing these results 
with accretion rates from Mendisutfa et al. who 
used a similar methodology, and from ICarcia Lopez et al.l 
l|2006h , who based their accretion rates on Br7 emission, 
provides literature values for 15 stars in our sample. Wor- 
ryingly, the literature accretion rates differ for five of these 
stars by between 0.5 and 1 dex. Whether this apparent dis- 
agreement is due to intrinsic stellar variability, differences 
in methodology, or sirn ply large uncertainties is not clear 
iMendigutfa et cil]|201ll ). We compare these accretion rates 
with our peculiarity index [P] (defined in Sect. [5|, stellar 
mass, and fractional pre-main sequence age in Fig. S) In 
cases where multiple accretion rates are available the aver- 
age value was taken. V380 Ori was left out of the comparison 
due to its binarity. We see no correlation between accretion 
rate and chemical peculiarity or mass. There is possibly a 
weak correlation with fractional pre-main sequence age, but 
due to the large uncertainties on both the accretion rates 
and the fractional ages it is not entirely clear. 

Our results provide further evidence against the hy- 
pothesis that A Boo stars are formed by radiative d i ffusio n 
with mass loss, as proposed by'Michaud & Charland) (|l986l ). 
In their calculations, iMi chaud & Charland (1986) find that 
it takes 10* or 10^ years to build up A Boo peculiarities by 
diffusion in the presence of mass loss. This is much greater 
than the ages of the A Boo-like stars in our sample, which 
are only a few x 10® years old. Thus the time scale for this 
mechanism is dramatically incompatible with our observa- 
tions. This is consistent with the generally accepted view 



that diffusion with m ass loss is insufficien t for producing A 
Boo peculiarities (e.g. ICharbonneau|[l993l ) . 



6.2 Magnetic stars 

The one possible Bp star in our sample is V380 Ori A. This 
star shows modest overabundances of Fe, Mn, and Ni of 
~0.6 dex, and Si is also enhanced relative to solar. C, Ne, 
and S are consistent with solar, as are O and Mg at 2a. This 
pattern is characteristic of Ap/Bp stars. Arguably, results 
are marginally consistent with the star simply having a high 
metallicity (~0.4 dex above solar). Ifowever, since Fe, Mn 
and Ni are more enhanced than 0.4 dex, and C, O, Ne, Mg 
and S are less enhanced than this, we consider it more likely 
that the star is weakly chemically peculiar. If the star is 
a Bp star, the star has begun developing Ap/Bp chemical 
peculiarities, but they have not yet reached the degree seen 
in most main sequence Ap and Bp stars. 

Atomic diffusion producin g chemical peculia rities was 
modelled in V380 Ori A by IVick et al.l (|201ll ). using a 
model that includes stellar evolution, as well as mass loss 
and atomic d iffusion. They used t he stellar parameters de- 
termined by lAlecian et all (|2009l ). who find a somewhat 
lower mass than we do (mostly due to assuming Ry ~ 3 
rather than 5, thus using less extinction). Nevertheless, the 
abundance s they derive should be roughly representative. 
IVick et al.l (|201ll ) find Mn, Fe and Ni to be overabundant 
by ~0.5 dex, which agrees nicely with our results. They find 
roughly solar abundances for C, O and S, which agrees with 
our results. The overabundance of P that they find is also 
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Figure 3. Trends in stellar parameters with the peculiarity index. The index is the average of Cr, Fe and Ni abundances each relative 
to solar, minus the average of the C, N and O abundances relative to solar. Plots are included for effective temperature, mass, age, and 
fractional pre- main sequence age (r). 



roughly consistent with our results. However, the under- 
abundant Mg, and solar Si abundance that they find are 
inconsistent with our results. Thus, while this model is im- 
perfect, we find it does provide a reasonable approximation 
of the majority of the observed peculiarities in V380 Ori A. 

Interestingly, HD 190073 and HD 101412 are both mag- 
netic but do not show Ap/Bp chemical peculiarities. HD 
190073 is chemically normal, while HD 101412 shows A Boo 
type chemical peculiarities. Thus, unlike on the main se- 
quence, not all pre-main sequence magnetic A and B stars 
are chemically peculiar. Atomic diffusion theory suggests 
that enough time must elapse after accretion and deep sur- 
face convection have halted for significant peculiarities to de- 
velop. The dissipation of a large convective envelope can be 
modelled easily enough, and in most models dissipates early 
in pre-main sequence life of A and B stars. However, deter- 
mining when accretion has dropped sufficiently for atomic 
diffusion to dominate is much harder. Indeed it may depend 
on the circumstellar environment of the particular star, and 
thus not be reliably predictable. 

The A Boo peculiarities of HD 101412 presumably sim- 
ply reflect the high incidence of A Boo peculiarities in HAeBe 
stars in general. Thus, we suspect that these peculiarities are 
not a consequence of the magnetic field of HD 101412. 

Modelling by IVick et al.l (|201ll ') suggests that HD 
190073 has not yet had time to develop Ap/Bp chemical 



peculiarities. In their models, a star of this mass requires 
~2 Myr to develop peculiarities by atomic diffusion, while 
we find HD 190073 is only 1.6to.6 old- HD 101412 has 
almost the same mass, and we find an age o f only 1.2^g * 
Myr, thus by the models of IVick et al.l (|201ll ) it is also too 
young to have developed Ap/Bp peculiarities. The younger 
age of HD 101412 is consistent with the picture of its A Boo 
type chemical peculiarities being due to selective accretion, 
since accretion is more likely to be ongoing earlier in a stars 
hfe. 

It is interesting that the two stars in this sample with 
the largest amount of emission in their optical spectra are 
V380 Ori A and HD 190073, both of which are magnetic. In 
particular, the emission mostly appears in lines that would 
be strongly in absorption for a purely photospheric spec- 
trum of the star. HD 101412, the third magnetic star, dis- 
plays a fairly modest amount of emission. Additionally, the 
measured accretion rates for V380 Ori and HD 190073 are 
both substantially higher than for other stars in our sam- 
ple (log Mace = -5.6 [Mayr-^] for V380 Ori, logM^cc = 
-5.0 [Mpyr"^] for HP 190073: [Poiiehew fc Brittainll201ll : 
iMendigutfa et ah! I2OIII ) . The impact of the secondary in 
V380 Ori does not appear to be considered in the mea- 
surement of the star's accretion rate, and thus it may be 
somewhat inaccurate. In classical T Tauri stars, the strong 
emission in their spectra is thought to result from magnet- 
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Figure 4. Trends in stellar parameters with literature accretion rates. Plots are included for our chemical peculiarity index, stellar mass, 
and fractional pre-main sequence age (r). 



ically channelled accretion flows falling onto the star. We 
speculate that the strong emission in V380 Ori A and HD 
190073 might be produced by a similar process. However, 
it is not clear whether this is the case, both because of the 
small sample size, and because the non-magnetic stars in the 
sample are biased towards lower emission stars (since large 
amounts of emission limit the abundance analysis). Mag- 
netic fields are not detected in individual emission lines for 
any of our stars, thus we cannot use those lines to assess the 
presence of magnetospheric accretion directly. 

The three stars in our sample with confirmed magnetic 
fields also have the lo west v sin i values (for more analysis see 
lAlecian et al.|[2"012bl ). In this respect these three magnetic 
HAeBe star are similar to Ap/Bp stars, which rotate slower 
than normal non-magnetic A and B stars. The low vsmi of 
the mag netic HAeBe stars H D 200775A (wsini = 26 ± 2 
kms-\ lAlecian et all l2008al) and of NGC 224 4 OI 201 
(vsmi = 23.5 ± 0.5 kms~\ lAlecian et al.ll2008bh . and the 
moderate v sin i of HD 72106 A (vsini = 41.0±0.6 kms~\ 
iFolsom et al.llioO^ ') support this observation. However, the 
very hi gh vsmi of NGC 6 611 W610 {vsmi = 190 ± 10 
km s~'^. lAlecian et al.|[2008bl ) indicates that not all magnetic 
HAeBe stars rotate slowly. Nevertheless, the large number 
of low vsini magnetic HAeBe stars suggests that mag- 
netic bra king may have occurred in many of these stars (see 
lAlecian e t al. 2012b). 

Several other very young magnetic A and B stars are 



known to have chemical peculiarities. iFolsom et all (|2008D 
examined the magnetic star HD 72106 A, and found Ap/Bp 
peculiarities. They find that the star has a mass of 2.4 Mq. 
While they were unable to firmly conclude that HD 72106 A 
is on the pre-main sequence, they do find that the system is 
between 6 and 13 Myr from the birthline, based on the H-R 
diagram position of the secondary in the binary system. The 
pattern of abundances seen in HD 72106 A is similar to that 
seen in V380 Ori A, but more extreme. Overabundances of 
iron peak elements and a weak overabundance of Si is seen 
in both stars. However, the strong He underabundance seen 
in HD 72106 A is not seen in V380 Ori A. The differences in 
abundances between the stars could be due to V380 Ori A 
being hotter and more massive HD 72106 A, or due to HD 
72106 A being more evolved. 

iBagnulo et al] (|2004D investigated the main sequence 
magnetic star NGC 2244-334, which is a member of the open 
cluster NGC 2244, and found Ap/Bp peculiarities. They 
found that the star has a mass of ~4 Mq and commented 
that, as a cluster member, the star is likely ~2 Myr old as 
measured from the ZAMS. This is a similar mass to V380 
Ori A, though NGC 2244-334 is somewhat older. The pecu- 
liarities they found are very similar to those of HD 72106 
A, with overabundances in iron peak elements of 1-2 dex, a 
somewhat weaker overabundance of Si, and a strongly un- 
derabundant He. This is again similar to what we see in 
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V380 Ori A, though much stronger, and with a He under- 
abu ndance. 

lAlecian et al.l (|2008bl ) studied the magnetic Herbig Be 
star NGC 6611 W601, and noted peculiarly strong He and Si 
lines in its spectrum. This suggest the star likely represents 
a pre-main sequence He-strong star, but a detailed abun- 
dance analysis has yet to be performed. While this pattern 
of peculiarities would be different from that of V380 Ori A, 
the star has a much higher temperature of 22500 rfc 2500 
K and a much higher mass of IO.2IJ;? Mm (|Alecian et al.l 
l2008bl ). For a main sequence magnetic star in this temper- 
ature and mass range, one would expect to see He-strong 
peculia rities, rather than Ap/Bp peculiarities. lAlecian et all 
l|2008bl ) estimate the age of NG C 6611 W601 to be less t han 
0.06 Myr from the birthline of IPalla fc Stahleil (|l993l . for 
an accretion rate of 10"* Mq yr~^), based on its H-R dia- 
gram positiom_WMet his is very young, there are no models 
from lVick et al.l (l201lh f o r a co mparable mass. We do note, 
however, that Vick et al] (|201ll ) find peculiarities developing 
faster for more massive models. 

Thus, while the peculiarities seen in V380 Ori A are 
weak, the pattern of peculiarities is consistent with those 
seen in other young intermediate mass magnetic stars. This 
is consistent with the hypothesis that the weak peculiarities 
we see in this star will develop into the stronger peculiarities 
seen in most main sequence Ap/Bp stars. 



6.3 Chemically normal stars 

We find no hints of any chemical peculiarities in the stars 
of our sample other than A Boo peculiarities and the weak 
Bp peculiarities in V380 Ori A. Thus we find no evidence 
for Am or HgMn stars on the pre-main sequence. Our sam- 
ple is not large enough to definitively exclude the existence 
of Am and HgMn stars at the same incidence as on the 
mai n sequence. However, if the roughly 10-20% incidence 
(e.g. ISmithlll996l l of these stars on the main sequence is ex- 
tended to the pre-main sequence, we would expect to have 
seen two to four such stars. Thus it is possible that these 
peculiarities do not form on the pre-main sequence, unlike 
for Ap/Bp stars. The atomic diffusion thought to give rise to 
these peculiarities may be disrupted by weak accretion more 
readily than the diffusion in the presence of a magnetic field 
that produces Ap/Bp stars. The presence of a magnetic field 
may stabilise a stellar atmosphere and increase diffusion ve- 
locities, potentially making diffusion in the presence of a 
magnetic field more robust (|Michaudlll970l :[ Michaud et al.l 
Il98ll ). Alternately, diffusion in non- magnetic stars may sim- 
ply require more time than the pre-main sequence lifetime 
of these stars to produce significant peculiarities. 

The large majority of young A and B stars analysed, 
both on the main sequence and on the pre-main sequence 
(e.g. lAcke fc Waelkensll2004 ). have been found to be chem- 
ically nor mal. A very precise s tudy of A and B star abun- 
dances by iFossati et al.l (|200 9l) investigated what 'normal' 
abundances are for A and B stars. They found that solar 
abundances are a good approximation, though there were 
small deviations for specific elements, which appeared to de- 
pend on the temperature of the star. St udies of nearby open 
clusters support thi s conclusion (e.g. Fossati et al.l I2OO8I : 



IVillanova et al.ll2009l:lFossati et alll201lf ) 

3r 



Acke fc Waelkena (120041 ) performed an abundance anal- 



ysis of 24 pre-main sequence and very young main sequence 
stars, based on equivalent widths. A comparison of individ- 
ual results for the 7 stars we have in common with their 
study is discussed in Sect. [X] Our results are generally con- 
sistent, although we are able to determine abundances for 
a wider range of elements, and in many cases provide more 
precise abundances. The one notable difference in our results 
is for HD 139614. For this star our mor e accurate TeS al- 
lows us to classify the star as A Boo, while lAcke fc WaelkensI 
(2004) could only comment that it a ppeared to be metal 
weak. In total I Acke fc Waelkend (|2004') found only one Her- 
big star that was also a clear A Boo star (HD 100546), and 
they noted that AB Aur shows hints of A Boo peculiar- 
ity. In part this difference in detection rates may due to 
our more accurate T^s and log g values, as the case of HD 
139614 demonstrates. Our higher precision, due to an im- 
proved methodology and higher S/N observations, is also 
a significant part of this difference. Finally, both our stud- 
ies deal with a fairly small number of stars, and thus the 
true incidence of A Boo peculiarities in HAeBe stars remains 
somewhat uncertain. 



6.4 Possible Veiling 

One concern when analysing spectra of T Tauri stars is veil- 
ing. Veiling is additional continuum emission, usually from 
an accretion shock (or other circumstellar material), which 
makes spectral lines appear weaker than they would in a 
purely photospheric sp ectrum. In HAeBe stars, v eiling can 
be seen m the UV fe.g. lDonehew fc Britta"inll201ll ). but usu- 
ally is considered to be unimportant in th e optical (e.g. 
'Bohm fc Catala'1993"; 'Herna ndez et al.ll2004l : ICowlev et al] 
2010; Donchc w fc Britt ain 201lh. HAeBe stars are intrinsi- 
cally much more luminous than T Tauri stars, thus much 
more emission from an accretion shock would be needed to 
have the same veiling effect. Measured accretion rates for 
HAeBe stars are usually not much larger than for T Tauri 
stars, thus it is unlikely that HAeBe stars would have much 
larger amounts of emission from accretion shocks. Conse- 
quently, we would not expect veiling to have much impact 
on HAeBe star spectra. 

If veiling is present in our observations it could make a 
star appear to be metal poor when it is not, thus it is impor- 
tant to check for veiling. If it is present in our spectra, veiling 
should have some wavelength dependence. We have com- 
pared abundances derived from different wavelength regions 
for individual stars. No clear correlation of abundance with 
wavelength was found, suggesting that veiling is unimpor- 
tant in these stars. Veiling should affect all lines in a wave- 
length region similarly. Thus, if it is present veiling should 
affect most elements similarly, since the lines for most ele- 
ments are widely distributed in our spectra. In particular, it 
could not produce underabundances in iron peak elements 
while keeping solar abundances of C and O that we see in 
many of the stars in our sample. Therefore veiling should not 
affect the relative abundances of C, N, and O to Fe, even if 
it were present at a significant level in our spectra. We con- 
clude it is unlikely that veiling has a significant impact on 
our derived abundances, and even if veiling is present it is 
not obvious how it could produce peculiar abundance pat- 
terns. 

Arguably, the pattern of A Boo underabundances we ob- 
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serve could be due to weak emission infilling in all the lines 
of some elements, however this is very unlikely. We have 
searched carefully for emission in lines, both by looking for 
unexpected variability and by looking for lines with shapes 
inconsistent with rotational broadening. Thus we have ex- 
cluded most if not all lines with emission. More importantly, 
lines of an element with lower excitation potentials are more 
likely to have stronger emission, thus emission would not be 
uniform across all lines. If emission infilling was occurring 
and otherwise undetected, we would observe it as an inabil- 
ity to simultaneously fit a wide range of lines of one element 
with one abundance. We do not find this, and consequently 
can safely conclude that the A Boo peculiarities we detect 
are real, and not simply due to emission infilling. 



7 CONCLUSIONS 

We have performed a detailed abundance analysis of 21 
stars, 20 of which are Herbig Ae and Be stars, 1 of which 
(tt Get) is not a Herbig star, but may still be on the pre- 
main sequence. Three of these stars have confirmed magnetic 
fields, while magnetic fields do not appear to be present at 
detectable levels in the remaining 18. 

We find 9 stars that are chemically normal (one of which 
is rather uncertain), 11 that show A Bootis chemical pecu- 
liarities, and one that is weakly Ap/Bp. This is a remarkably 
large fraction of stars that display A Boo type peculiarities. 
On the main sequence, onl y ~2% of stars display A Boo 
peculiarities (jPaunze We interpret this as evidence 

in favour of a selective accretion hypothesis for the devel- 
opment of A Boo peculiarities. Since Herbig stars are either 
accreting or have recently finished accreting, a process which 
depends on accretion should be seen more frequently in these 
stars. 

Of the three magnetic stars, HD 101412 displays A Boo 
peculiarities, HD 190073 is chemically normal, and V380 Ori 
is weakly Bp. Thus it appears that magnetic Herbig stars 
can display Ap/Bp type peculiarities, but not all magnetic 
Herbig stars do. This is in contrast to main sequence mag- 
netic A and l ate B stars, which a lways display Ap/Bp type 
peculiarities (|Auriere et all bOOTh . The presence of A Boo 
peculiarities in a magnetic star appears to simply reflect the 
high incidence of these peculiarities in all Herbig stars. 
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APPENDIX A: RESULTS FOR INDIVIDUAL 
STARS 

Non-magnetic stars 
Al TT Get (HD 17081) 

TT Get (HD 17081) i s a coo l B star that was studied in de- 
tail bv lFossati et al.l (| 20091 ). They performed a very precise, 
detailed abundance analysis, and considered the star to be 
a good example of a 'normal' cool B star, tt Get is not a 
Herbig Be star by the strict definition, since it has no strong 
emission in its optical s pectrum. However, the star does dis- 
play an infrared excess jMalfait et aL iggsl). as w ell as some 
unexplained line profile variations IjFossati et al.. 2009), and 
an H-R diagram position to the right of the main sequence. 
Thus the star is certainly very young, and possibly still on 
the pre-main sequence, despite the lack of clear emission. 
Based on the H-R diagram position of the star, it appears 
to be less than 0.5 Myr from the birth line , and almost cer- 
tainly is less than 1 Myr old. lMalfait et al] (|l998l ) suggested 
this star was a Vega-type star, a classification which is con- 
sistent with our observations. 

In order to check the accuracy of our methodology, we 
analysed an observation of tt Get an d com pared our re- 
sults with the results of iFossati et al. 



( 20091). Our result s 



are completely consistent with those of lFossati et al.l (|2009l ). 
for temperature, gravity, t;sini, and abundances. Our atmo- 
spheric parameters and abundances are compared with those 
of Fossati et al. (2009) in Table ful and the differences be- 
tween our abundances and those of Fossati et al.l l|2009l ) are 
plotted in Fig. IA1[ The excell ent agreement between our re- 
sults and lFossati et af] (|2009l ) demonstrates the accuracy of 
our methodology. 

TT Get was observed spectropolarimetricly by 
IWade et all (|2007l '). who did not detect a magnetic 
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Figure Al. The di fference between our final abundances and 
those o f iFossati etH ( 2009) for vr Get. A good agreement between 
the two sets of abundances is found. 



Table Al. Our results compared 
||2009| ') for TT Get. 



with those of iFossati et al.l 





Fossati et al. (2009) 


This work 


Toff (K) 


12800 ± 200 


12900 ± 400 


log 9 


3.75 ±0.1 


3.8 ±0.2 


V sini (km s~-^) 


20.2 ±0.9 


20.9 ± 1.2 


C (kms-i) 


1.0 ±0.5 


1.7 ± 1.0 


He 


-0.97 ±0.04 


-0.91 ±0.14 


C 


-3.58 ±0.07 


^ -3.4 


N 


-4.03 ±0.13 


-4.03 ±0.15 


O 


-3.06 ±0.14 


-3.16 ±0.13 


Ne 


-3.66 ±0.09 


-3.76 ±0.28 


Na 


-5.23 ±0.07 


-5.3 ±0.1 


Mg 


-4.47 ±0.16 


-4.43 ±0.15 


Al 


-5.73 ±0.27 


-5.81 ±0.17 


Si 


-4.41 ±0.2 


-4.55 ±0.14 


P 


-6.38 ±0.19 


-6.43 ±0.14 


S 


-4.78 ±0.16 


-4.99 ±0.08 


Ar 


-5.24 ±0.19 


-5.6 ±0.3 


Ca 


-5.77 


-5.7 ±0.2 


Ti 


-7.42 ± 0.08 


-7.37 ±0.15 


Or 


-6.41 ±0.1 


-6.57 ±0.16 


Mn 


-6.5 ±0.09 


-6.6 ±0.3 


Fe 


-4.58 ±0.14 


-4.70 ± 0.08 


Ni 


-5.76 ±0.19 


-5.85 ±0.06 



field. lAlecian et al.l l|2012ah also found no magnetic field in 
TT Get, and thus we consider the star to be non-magnetic. 

A2 HD 31293 (AB Aur) 

We derive T^a = 9800 ± 700 K and log g = 3.9 ± 0.3 
for HD 31293 (AB A ur). This T^s is consistent with 
lAcke fc Waelkena l|2004 ). but our log g is inconsistent with 
their (remarkably hi gh) value of 5.0. The T^s and log g axe 
both consistent with iBohm fc Cataial (|l993l ). We find solar 
He and O abundances, C and N are la above solar, while 
the iron peak elements are consistently la below solar. The 



systematic underabundance of iron peak elements by ~0.3 
dex suggests the star may be a weak A Boo star. This is 
strengthened somewhat by the observation that the iron 
peak elements are ~0.5 dex below the average C, N, and O 
abundance. Based on these underabundances, we conclude 
that the star most likely is a weak A Boo star, but it is a 
mar ginal case. 

lAcke fc Waelkenj l|2004l ) derive abundances for O, N, 
and Si that are consistent with ours, but their Fe abundance 
of 1 dex below solar (based on 2 Fe I lines) is inconsistent 
with our value. This d i fferen ce may be due to the higher 
log g lAcke fc Waelkeni ()2004l ) used. Alternately, it is possi- 
ble their measurement was impacted by the large amount of 
emission in the spectrum of HD 31293. Both Fe abundances 
are consistent with a A Boo classification of the star, but 
we prefer our larger abundance since it is based on larger 
number and a wider range of lines, and it is consistent with 
the neighbouring i r on pea k element abundances we derive. 

[Alecian ct al.' ('2012a') find no magnetic field in HD 
31293. Wade ct al. (2007) also find no magnetic field in this 
star, thus we consider it to be non-magnetic. 



A3 HD 31648 

For HD 31648 we derive an T^a of 8800±190 K and a log g of 
4.1 ±0.2, which is consistent with Ackc fc Waclkcns (2004). 
We find most abundances within la of solar, and all abun- 
dances except Mg and Ti within 2a. The abundances are 
systematically slightly above solar and may refiect a slightly 
enhanced metallicity, or a slightly overestimated Tcff. 

Our Ca and Fe abundances are consistent with 
lAcke fc Waelken^ (|2004l ). our C abundance is consistent at 
2a, as is our Y abundance if we assume the uncertainty on 
their value is comparabl e to ours. Our Si an d Mg a bundances 
are not consistent with lAcke fc Waelkenj (|2004l ). although 
their abundances are based on only 1 line and thus are likely 
uncertain. The large vsmi of this star causes significant 
blending in its spectrum. This co uld explain the discrepan - 
cies between our results, since the lAcke fc Waelkend (|2004l ) 
abundances are based on equivalent width measurements, 
which can easily be influenced by unaccounted-for blending. 

HD 31648 was not found to have a magnetic held 
by 'Alecian et al.' (' 2012ah . No magnetic fleld w as found by 
IWadc ct al. (2007) either. iHubrig et all (|2006l ) claimed a 
weak detection of a magnetic fleld in their observation of HD 
31648 but, in the r e-reduction of the observation made by 
iHubrig et all (|2007|) the field is below their detection thresh- 
old. Surprisingly, Hubrig et al.l (|2007l ) report a signal the H/3, 
H7, H(5, and Ca H & K lines of their V/I spectrum, which 
they interpret as a "circumstellar" magnetic field. This re- 
sult has yet to be independently confirmed. A careful exam- 
ination of the H/3, H7, H(5, and Ca H fc K lines in ou r spec - 
tra shows no polarisation signature. Bagnulo e t al.l (|201ll ) 
re-reduced a wide range of spectropolarimctric observations 
from the FO RSl instrument at the VLT, including the ob- 
servations of lHubrig et al.l (|2006l ) . They find no evidence for 
a magnetic field in HD 31648. Considering the lack of evi- 
dence for a photospheric magnetic field, and the absence of 
a confirmation of a "circumstellar" magnetic field, we con- 
clude that HD 31648 is non-magnetic. 
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A4 HD 36112 

For HD 36112 we find T^a = 8190 ± 150 K and log g = 
4.1 ± 0.4, both of which are so mewhat greater than th e 
Teff = 7750 and log g = 3.5 of lAcke fc Waelkend (|2004h . 
Our uncertainty on TeS may be overly optimistic, though it 
does refiect the scatter in Tcs derived from different spec- 
tral windows. The uncertainties on the calibration of Tefr 
become important at this level, and thus there may be some 
systematic uncertainty not included in our error bar. Most 
of our derived abundances are within la of solar, and all 
are within 2a. There are no systematic enhancements or de- 
pletions, thus we conclude HD 36112 is chemically normal, 
with solar abundances. 

Our abundances are consistent with those of 
lAcke fc Waelkend (|2004l ) for C, Si, S, Ca, Sc, Cr, Mn, 
Fe, and Ni. Our abundance for Ti and Y differ by ~2cr. 
These differences may be a consequence of our hotter 
temp erature, or they may be due to the small number of 
lines lAcke fc Waelkend l|2004l ) used for determining those 
abunda nces. 

Alecian et al.l ()2012al ) do not detect a magnetic field in 



this star. In their FORSl observations Wade et al 



do not detect a magnetic field, nor do lBagnulo et al 



2007) 



2011 



in their re-reduction of the FORSl observations. 



A5 HD 68695 

For HD 68695 we derive T^s = 9000 ± 300 K and log g = 
4.3 ± 0.3. We find C, N, O, and S abundances that are ~lcr 
above solar, while iron peak abundances are 2 to 3a below 
solar (~0.7 dex below solar), except for the rather uncertain 
Mn abundance. We conclude that this star has clear A Boo 
peculiarities. 

No magnetic fiel d was detected in HD 68 695 by either 
lAlecian etall (12012^ 1 or bv lWade et al.l l|2007h . thus we con- 
sider the star to be non-magnetic. 



A6 HD 139614 



HD 139614 was analysed bv lAcke fc Waelkend l|2004h using 
an Tcff of 8000 K. They found an overall deficiency of metals, 
but no clear evidence of the selective depletion characteristic 
of A Boo stars. We find HD 139614 to have an T^a of 7600 ± 
300 K, and a log g of 3.9 ± 0.3, which is consistent with the 
res ults of Teff = 7400 ±200 K and log g = 4.0 ±0.4 reported 
by Guimaraes et al.l (l2006f ). and marginally consistent with 
lAcke fc Waelkend (|200^ fat 1.3ff in Teg and 2a in log g). 
We consider our para meters (and Guimaraes et al.ll2006 ) to 
be more accurate, since^Ac ke fc Waelk ens (2004) values were 
based on the spectral classification of Malfait ct al. (1998), 
which could be infiuenced by chemical peculiarities in the 
star. In fact, for our parameters we find a clear A Boo pattern 
of abundances, with solar C, N, and O, and underabundant 
iron peak elements. 

No magnetic fi el d was detected in HD 139614 by 
lAlecian et all (|2012al ). iHubrig et al.l lj2004h claimed a detec- 
tio n of a magnetic fie ld in this star, which was not supported 
bvlWade et all l|2005l ). Further observations bv lHubrig et al.l 
(2009") did not detect a magnetic field in the star, though 
tliey comment on a possible signal in the Ca H and K lines . 
In their re-reduction of FORS data, iBagnuto et al.l (120111 ) 



find no evidence for a magnetic field i n HD 139614, and co n- 
clude that the reported detection bv lHubrig et al.l (|2004l ) is 
likely spurious. Consequently, we consider HD 139614 not 
to have a significant magnetic field. 



A7 HD 141569 

We find Toff = 9800 ± 500 K and log g = 4.2 ± 0.4 for HD 
141569, which is consis tent with the parameters found by 
iGuimariies et all l|2006i ). We find A Boo peculiarities in this 
star, with solar He, C and N, possibly overabundant O, and 
underabundant iron peak elements. 

Alecian ct al. (2012 j) do not detect a magnetic field in 
the star. Wade et al.l (j2007l ) also find no magnetic field, thus 
we conclude HD 141569 does not have a significant magnetic 
field. 



A8 HD 142666 

For HD 142666 we find T^b = 75 00 ± 200 K and l og g = 
3.9 ± 0.3, which is consistent with IGuimaraes et al.l (|2006D . 
We find weak A Boo peculiarities in this star. C, O, Na 
and S have solar abundance, while the iron peak elements 
are consistently ~0.5 dex below solar, except for the rather 
unc ertain Co ab undance . 

IWade et al. (2007) observed the star and found no mag- 
netic field. Alecian et al.l (|2012al ) find no magnetic field in 
the star either. 



A9 HD 144432 

For HD 144432 we find T^s = 74 00 ± 200 K and l og g = 
3.9 ± 0.3, which is consistent with IGuimaraes et al.l (|2006h . 
The abundances of all elements we derive are within 2a of 
solar, and most are within la. There are no clear trends in 
abundances of different elements, thus we consider the star 
che mically normal , with n early solar metallicity. 

lAlecian et al] (|2012al ') observed HD 144432 twice and 



found no magnetic field IWade et all tOOlh also found no 
magnetic field in the star. iHubrig et aTr(|2007l ) claim a 
weak magn etic field dete c tion in thi s star in one obser- 
vation, but IHubrig et all (120041 ) and IHubrig et all (|2009l ) 
did not detec t a s ignificant magnetic field in the star. 
iBagnulo et al] (|201ll ) re-reduced the one observation of HD 
144432 with a magne t ic fiel d detection from IHubrig et alJ 
(|2007l ). [Bagnulo et all (|201ll ') do find a magnetic field in the 
observation at a 4.9a level {{B^) = -108 ± 22 G). How- 
ever, they caution that the impact of instrumental effects 
on the detection of weak magnetic fields is not well under- 
stood for the FORSl instrument. They conclude that for 
weak magnetic fields, detections at less than 6a should not 
be considered definite, and require supporting observations. 
Thus, while we cannot completely exclude the possibility of 
a magnetic field in HD 144432, based on the large majority 
of the evidence (five out of six observations) we conclude 
that HD 144432 likely has no significant magnetic field. 



AlO HD 163296 

We determine an Tcff of 9 200 ±300 K and a l og g o f 4.2±0.3, 
which is consistent with IGuimaraes et al] (|2006l ). We find 
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the abundances of most elements to be within uncertainty of 
solar, though Mg and Cr are both overabundant at slightly 
more than 2a. With largely solar abundance and no clear 
selective depletion or enhancement, we conclude that the 
star is chemically normal . 

lAlecian et al.l ('2012a'') do not find a magnetic field in HD 
163296. iHubrig et al., (.2009 ) also do not detect a magnetic 
field in the star, thus we conclude the star is non-magnetic 



All HD 169142 

For HD 16914 2 we find T^s ^ 7500 ± 200 K, which is 
consistent with iGuimaraes et al.l l(2006l), and log g = 4.3 ± 
0.2, which is 2a greater than IGuimaraes et al. 1 (|2006l ) who 
find log q — 3.7 ± 0.1 . If we were to adopt the log g of 
IGuimaraes et all (|2006l ) . it would decrease our abundances 
by between ^0.05 and ~0.2 dex depending on the element, 
which is generally within our uncertainties. We find a clear 
A Boo pattern of abundances in this star. C and O are so- 
lar, and S is only slightly underabundant, while iron peak 
elements are underabundant by between 0.5 and 1 dex. In- 
terestingly Ba and Nd are roughly solar, though Eu may be 

und erabundant . 

lAlecian eTZI tOUsj) do not detect a magnetic field in 
HD 169142, nor do lHubrig et"all (|2009l ). 



A12 HD 176386 

We derive an T^a of 11000 ± 400 K and a log g of 4.1 ± 
0.3 for HD 176386. The abundances we derive for the star 
are consistent with solar, with the possible exception of Ca. 
Thus we conclude that the star is chemically normal. 

No magnetic fi e ld was found for HD 176386 by 
lAlecian et al.l l|2012al ') . iHubrig et al.l (|2009|) claim a margina l 
detection at 3a in one observation. iBagnulo et"al] (|201ll ) 
re-reduced this observation and found no magnetic field. 
Thus, based on the majority of the evidence, we consider 
HD 176386 to be a non-magnetic star. 



A13 HD 179218 

For HD 179218 we derived Toff = 9640 ± 2 50 K and log g = 
3.9± 0.2, both of which are consistent with IGuimaraes et al.l 
ll2006il . We find a solar He abundance, abundances of C, 
O, and Na that are marginally enhanced relative to solar, 
and N and S abundances that are uncertain, but apparently 
enhanced relative to solar. Iron peak elements are depleted 
relative to solar by ~0.5 dex. We conclude that HD 179218 
displays A Boo peculiarities. The overabundances of C, N, O, 
and S may reflect an above solar intrinsic metallicity in the 
star, which the A Boo peculiarities are superimposed upon. 
However, this is not certain, as He appears to have a nearly 
sola r abundance. 

lAlecian et al.l l|2012ah find no magnetic field in HD 

179218, whi ch is supported by the non-detection of 

IHubrig et all (|2009l ). 



A14 HD 244604 

We found Tofi = 8700 ± 220 K and log g = 4.0 ± 0.2 for 
HD 244604, which is consistent with the parameters from 



lAcke &: Waelkenj (120041 ). We find most elements have so- 
lar abundances or are enhanced by around la. There is no 
evidence for selective depletion or enha ncement. We con- 
cludc that the star is chemically normal. lAcke fc Waelkeii3 
(2004) were only able to determine an iron abundance for 
HD 244604, though it is consistent with our abundance if 
we assume their uncertainty is similar to ours. 

No magnetic field in HP 2446 04 is found by 
lAlecian et al] (|2012al ). IWade et al.l (|2007l ) also find no mag- 
netic field in the star, thus we consider the star to be non- 
magnetic. 



A15 HD 245185 

For HD 245185 we derive T^g = 9500 ± 750 K and l og g = 
4.0 ± 0.4, which is consistent with lAcke fc WaelkensI (|2004 ). 
We find He, C, N, and O abundances that are consistent 
with solar, and iron peak abundances that are ~0.8 dex 
below solar. While the uncertainties are relatively large for 
this star, the stro ng iron peak undcrabundances indicate A 
Boo peculiarities. lAcke fc Waclkcns. (.2004 ) were unable to 
determine any abundances for this star due to a low S/N in 
thei r observation and the high wsini of the star. 

lAlecian ^TZI (l2012a^ do not detect a magnetic field in 
the star, nor do lWade et al.l (|2007l ). Consequently, we con- 
sider the star to be non-magnetic. 



A16 HD 278937 (IP Per) 

For HD 278937 (IP Per) we find T^s = 8000 ± 250 K 
a nd log g = 4.1 ± 0.2, which is consistent with the T^b 
of iMiroshnichenko et al.l (|200ll ). though they derive a log g 
of 4.4 ±0.1 which is l a above our value . The star is a 
known 5 Scuti pulsator (|Ripepi et al.l [20061 ). We find solar 
abundances for C, N, O, and S, while iron peak elements 
as well as Na, Mg, and Si are between 0.5 and 0.7 dex 
underabundant. This star sh ows clear A Boo peculiarities. 
IMiroshnichenko et al.1 (|200ll ) found that HD 278937 had a 
metallicity below solar {[M/H] = —0.4), based on a fit to 
the star's metallic line spectrum. This value is most likely 
a result of the low iron peak abundances in the star, and 
thus consistent with our abundances. The solar abundances 
of C, N, O and S were likely overlooked by t he sca led so- 
lar abundance model of lMiroshnichenko et al.l ((SooJ) due to 
the relative scarcity of lines of those elements. Therefore, 
the star is not a low metallicity star, but rather appears to 
have peculiar photospheric abundances. 

No magnetic fi eld is found in HD 278937 by 
lAlecian et atl (l2012al ). or bv lWade et ai] (|2007l ). 



A17 T Ori 

We find Tcff = 8500 ± 300 K and log p = 4.2 ± 0.3 for T 
Ori. The abundances for He, C, N, and O are all consistent 
with solar, while the S abundance is almost 3a above solar. 
The S abundance is based only on two weak lines at 6749 
and 6757 A, and thus may be somewhat more uncertain 
than the error bar suggests. The iron peak abundances are 
clearly ~0.5 dex below solar. We conclude that T Ori has 
clear A Boo peculiarities 
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Neither lAlecian et al] l|2012al ') nor IWade et al] (|2007l ') 
detect a magnetic field in T Ori. 



Magnetic stars 



A18 HD 101412 



HD 101412 was observed by I Wade et al] (|2005l ) and reanal- 
ysed bv lWade et al. I (|2007l ) who found a longitudinal mag- 
netic field of 500 ± 100 G. The presence o f a magneti c field 
was recently confirmed by ^Hubrig et al.' (20091, 1201 ll ). and 
Alecian et al. (in prep.). I Cowley ct al. (2010) derived chem- 
ical abundances for HD 101412, and surprisingly discovered 
it had A Bootis peculiarities. In our abundance analysis we 
find underabundances of iron peak elements, and solar C, 
O, and He, with va lues that are consistent with those of 
ICowlev et al] l|2010i ). Thus we confirm that HD 101412 has 
A Boo peculiarities. While in main sequence A and B stars 
magnetic fields are associated with Ap and Bp chemical pe- 
culiarities, that is not the case for HD 101412. 

We confirm the presence of an "a nomalous saturation " 
in the lines of HD 101412 noted bv ICowlev et all l|2010l ). 
Even with a microturbulence of kms~^, when we fit weaker 
lines, stronger lines in the synthetic spectrum are often too 
deep to fit the observation. This problem is not consistent 
across all lines, and often appears to be fairly minor. This 
anomalous saturation does not appear to be common among 
HAeBe stars, as this is the only object for which we have 
seen this problem. A possible explanation for the discrep- 
ancy between the observations and the synthetic spectra is 
the presence of chemical stratification in the stellar atmo- 
sphere. Since stronger lines are often formed higher in the 
stellar atmosphere, if the abundance of some elements was 
lower higher in the atmosphere that would explain the ap- 
parent weakness of many stronger lines. The possibility of 
stratification will be investigated in detail in a future paper. 
The presence of veiling in the spectrum of this star would 
produce a similar discrepancy to what we see, however veil- 
ing would affect all lines of the same strength similarly, un- 
like the non-uniform trend we see. Thus we conclude veiling 
is un likely to be th e source of this discrepancy, in agreement 
with ICowlev et all (|2010l ). The discrepancy is not likely to 
be a direct effect of the star's magnetic field on line forma- 
tion, which usually desaturates lines, and we do not find the 
problem in the other magnetic stars of our sample. 

Placing HD 101412 on the H-R diagram, we find it has 
an age of 1.21q7 Myr and a mass of 3.0 ± 0.3 Mq, both 
of which are consistent with IWade et al.1 (|2007l ). Since the 
star has no Hipparcos pa rallax, we adopted the distance of 
500-700 pc proposed bv IVieira et al.l ()2003h . however this 
is somewhat uncertain. Even with this uncertainty, the star 
falls well before the beginning of the main sequence on the 
H-R diagram, thus the star is almost certainly a pre-main 
sequence object. 

We find similar A Boo type peculiarities in many of the 
non-magnetic HAeBe stars in our sample, thus this feature 
is not unique to HD 101412, or to magnetic HAeBe stars. 
Indeed, it appears that the peculiarities in HD 101412 occur 
in a significant fraction of all HAeBe stars. 



A19 HD 190073 

HD 190073 has a magnetic fi e ld, d iscovered and repeat- 
edly detected by ICatala et al.l (|2007h , with a longitudinal 
strength of 74 ± 10 G that was constant over the course 
of their observations. In our analysis we find that the star 
is chemically normal, as shown in Fig. IA5I Apparent over 
abundances relative to solar in N and Na may be due to 
non-LTE effects, or weak undiagnosed blends, since for both 
elements only a couple of very weak lines were usable for our 
ana lysis. 

lAcke fc WaelkensI (|2004l ') analysed HD 190073 and 
adopted T^a — 9250 K and log g — 3.5, which is consistent 
with our values of Tcff = 9230 ±260 K and log g = 3.7 ±0.3. 
Based on equivalent widths, they found abundances for C, 
N, O, Mg, Si, Ca, Sc, Ti, V, Cr, Fe, Ni Y, Zr, and Ba, all of 
which are consistent with ours. For Al they found an abun- 
dance -0.45 dex below solar, based on only one line, which 
is inconsistent with our value. Our abundance is based on 
two lines, thus we consider it to be the more accurate value. 
Our result leads to a nearly solar abundance, consistent with 
other elements in this star. 

We find that the star is quite you ng (1.6 ± 0.6 Myr) , 
which agrees with the age determined bv lCatala et al.l (|2007| ) 
of 1.2 ±0.6 Myr. Thus it is possible that the star has not yet 
had time to build up significant chemical peculiarities on its 
surfac^^ Th is hyp othesis is supported by the modelling of 
IVick et al] (|201ll l. who found that it takes at least ~2 Myr 
for a 2.9M0 star to develop significant chemical peculiarities 
by radiative diffusion. 



A20 V380 Ori A & B 

V380 Ori is a triple star s ystem, with two com ponents de- 
tectable spectroscopically l| Alecian et alj|2009l ). The more 
massive spectroscopic component (the primary) shows no 
detectable change in radial velocity, while the fainter, less 
massive (secondary) co mponent's radial ve locity varies with 
a period of ~100 days (| Alecian et al.ll2009D . The third com- 
ponent in the sys tem was detected by speckle interferometry 
in the infra red (jLeinert et al 

][l993), 

but is not detectable 
in the visible spectrum of the system ( jAlecian et al. 200Sj[). 
The p rimary has a magnetic field, detected by iWade et alj 
(|2005l ), while the secondary has no detectable magnetic field. 
Th e magnetic fie l d of th e primary was repeatedly observed 
bv lAlecian et ID 11200911 . who modelled the geometry as a 
dipole, with a strength of 2120 ± 150 G, and an obliquity of 
66 ±5°. 

Modelling the spectrum of the V380 Ori system was 
done somewhat differently than for the other stars in this 
study. Spectral disentangling was inapplicable due to the 
large amounts of variable emission in the spectrum, and the 
lack of radial velocity variation in the primary. Instead we fit 
a composite synthetic spectrum directly to the observation. 
The composite spectrum was constructed by adding the syn- 
thetic spectra of both stars, in absolute flux, weighted by the 
squared ratio of radii of the stars (the T^g dependence is in- 
cluded in the absolute fluxes of the synthetic spectra). The 
combined spectrum was then normalised by the sum of the 
continuum fluxes, also weighted by the ratio of radii squared. 
Using this composite spectrum, the model parameters for 
one star were simultaneously flt by minimisation, while 
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the other star's parameters were held fixed. The ratio of 
radii was simuhaneously fit as weU, but the radial velocities 
of the two components were determined before the rest of 
the fitting process. The combination of parameters for both 
stars from one spectral window was determined iteratively, 
alternating between the two stars, and fitting the parame- 
ters of one while holding the para meters of the o t her fi xed. 
Initial parameters were taken from lAlecian et all (|2009l '). 

We find chemical peculiarities in the primary that are 
consistent with Ap/Bp peculiarities. The peculiarities are 
weak, but detectable. C, S, and Ne are within uncertainty 
of solar, while Fe, Ni, and Mn are overabundant by more 
than 3(T. Si is also overabundant by more than 4(7 in the 
primary, and He may be overabundant at 2.2(t. The abun- 
dances of the secondary are uncertain, but consistent with 
solar. Most elements in the secondary are within Icr of so- 
lar, and all elements are within 2a (except for Ba). The large 
uncertainty in the secondary stems mostly from the uncer- 
tainty in the ratio of radii (and hence ratio of luminosities) 
for the system, though uncertainties in effective tempera- 
ture also play a role. While the average abundance in the 
secondary appears to be slightly above solar, it is not clear 
whether this reflects an intrinsically higher metallicity in the 
system, or if the luminosity ratio has simply been slightly 
overestimated. The abundances derived for the secondary 
are highly sensitive to the adopted ratio of radii, however 
the abundances for the primary are not, since the primary 
dominates the luminosity of the system. 

We find a ratio of radii for the V380 Ori system of 
Ra/Rb = 0.7 ± 0.1. The relatively large radius of the sec- 
ondary is a consequence of the star being less massive, and 
thus being slower to contract towards the main sequence 
than the more massive primary. The H-R diagram positions 
of the two stars are consistent with a single isochrone. 

The appearance of chemical peculiarities in t he primary 
of V3 80 Ori is consistent with the modelling bv IVick et al.l 
who discuss this system in some detail. Since the 
primary is more evolved than HD 190073, it could have de- 
veloped chemical peculiarities while HD 190073 has not yet. 
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Figure A2. Final abundances relative to solar for the stars in this study. Solar abundances are taken from lGrevesse et al.l l l2005h . Points 
marked with an arrow represent upper limits only. 
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Figure A3. Final abundances relative to solar for the stars in this study, as in Fig. IA2| 
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Figure A4. Final abundances relative to solar for the stars in this study, as in Fig. IA2I 
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Figure A5. Final abundances relative to solar for the stars in this study, as in Fig. IA2| 
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Figure A6. Comparison of the observed spectrum (jagged line) to the best fit synthetic spectrum (smooth line) for tt Get (HD 17081). 
Two independent wavelength regions are presented. Lines have been labelled by their major contributing species. Gaps in the synthetic 
spectra indicate regions that were not fit. 
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Figure A7. Gomparison of the observed spectrum to the best fit synthetic spectrum for HD 31293, as in Fig. IA6l 
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Figure A14. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 144432, as in Fig. IA6| 
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Figure A16. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 169142, as in Fig. IA6| 
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Figure A17. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 176386, as in Fig. IA6| 
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Figure A18. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 179218, as in Fig. IA6| 
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Figure A19. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 244604, as in Fig. IA6| 
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Figure A24. Comparison of the observed spectrum to the best fit synthetic spectrum for HD 190073, as in Fig. IA6| 
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Figure A25. Comparison of the observed spectrum (jagged line) to the best fit synthetic spectrum of the combined V380 Ori A and B 
system (smooth line). The upper two rows of line labels (green) indicate lines of the primary, the lower two rows (blue) indicate lines of 
the secondary. Gaps in the synthetic spectrum indicate regions that were excluded from the fit. 
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